Nourdin-Peccati analysis on Wiener and Wiener- Poisson space for 

general distributions 

Richard Eden Juan Viquez 

Department of Mathematics, Purdue University 
150 N. University St., West Lafayette, IN 47907-2067, USA 



Abstract 

Given a reference random variable, we study the solution of its Stein equation and obtain universal 
bounds on its first and second derivatives. We then extend the analysis of Nourdin and Peccati by 
bounding the Fortet-Mourier and Wasserstein distances from more general random variables such as 
members of the Exponential and Pearson families. Using these results, we obtain non-central limit 
theorems, generalizing the ideas applied to their analysis of convergence to Normal random variables. 
We do these in both Wiener space and the more general Wiener-Poisson space. In the former, we study 
conditions for convergence under several particular cases and characterize when two random variables 
have the same distribution. As an example, we apply this tool to bilinear functionals of Gaussian 
subordinated fields where the underlying process is a fractional Brownian motion with Hurst parameter 
bigger than 1/2. In the latter space we give sufficient conditions for a sequence of multiple (Wiener- 
Poisson) integrals to converge to a Normal random variable. 

1 Introduction 

Recent years have seen exciting research on combining Stein's method with Malliavin calculus in proving 
central and non-central limit theorems. The delicate combination of these tools can be attributed to Nourdin 
and Peccati who intertwined an integration by parts formula from Malliavin calculus with an ordinary 
differential equation called a Stein equation. Much work has been done to compare Normal or Gamma 
random variables (r.v.'s) with another r.v. (having unknown distribution). See |13j . |16j . |20j . [21] for results 
on the convergence of multiple (Wiener) integrals to a standard Normal or Gamma law. [3] and [37] discuss 
Cramer's theorem for Normal and Gamma distributions applied to multiple integrals. |29| gives probability 
tail bounds in terms of the Normal probability tail, with [5] applying the same techniques to give tail bounds 
in terms of the probability tail of other r.v.'s (e.g. Pearson distributions). 

In [T4], Nourdin and Peccati found a clever link between Stein's method and Malliavin calculus. This 
was used to derive the Nourdin-Peccati upper bound (NP bound) on the Wasserstein, Total Variation, 
Fortet-Mourier and Kolmogorov distances of a generic r.v. from a Normal r.v., and lay the groundwork 
for comparisons to a more general r.v. (with such results leading to non-central limit theorems). These 
authors and Reinert (see [17) ) applied this NP bound to obtain a second order Poincare-type inequality 
useful in proving central limit theorems (CLTs) in Wiener space. Specifically, they proved CLTs for linear 
functionals of Gaussian subordinated fields. Particular instances are when the subordinated process is 
fractional Brownian motion (fBm) or the solution to the Ornstein-Uhlenbeck (0-U) SDE driven by fBm. 
They also characterized convergence in distribution to a Normal r.v. for multiple stochastic integrals. 
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Later in (32] these ideas were applied to prove the NP bound in Poisson space (pure jump processes), 
which was used to obtain Berry-Esseen bounds for arbitrary tensor powers of 0-U kernels. Keeping in line 
with attempts to extend these results as far as possible, [3D] proved an NP bound in Wiener-Poisson space. 
The author applied similar ideas found in [17] to derive a second order Poincare-type inequality and use it 
to prove CLTs for a continuous average of a product of two 0-U processes (one in Wiener space and the 
other in Poisson space) which lives in the second chaos of Wiener-Poisson space. Also, it was proved that 
under mild conditions, the small jumps part of a functional in the first Poisson chaos is approximately equal 
in law to a functional in the first Wiener chaos with the same kernel (useful when simulating a fractional 
Levy process as a process with finitely many jumps plus a fBm). All these results show the importance of 
this NP bound and the potential it has as an effective tool in proving non-central limit theorems, CLTs and 
characterizations. 

Let Z be absolutely continuous with respect to Lebesgue measure. For our purposes, we can think of Z 
as a "well-behaved" r.v. (e.g. it lives in Wiener space and we know its density). Typical instances arc when 
Z is Normal, Gamma, or another member of the Pearson family of distributions. X is another r.v. whose 
properties are not as easy to determine as with Z, our "target" r.v.. We may have a hunch that X has the 
same distribution as Z^ or in the case of sequences, a belief that {X^^ converges in law to the distribution 
of Z. We thus want to compare X with Z. How different are the laws of X and Z for instance (and we need 
to make precise the sense in which they are different)? What conditions will ensure that X has the same 
law as Z1 For a sequence {A^n}, what sufficient conditions ensure convergence to Z in distribution? In this 
regard, we wish to measure the distance between (the laws of) X and Z hy a metric d-u which induces a 
topology that is equal to or stronger than the topology of convergence in distribution: if d-u (A„, Z) — > 0, 
then Xn — > Z in distribution. 

The motivation for this paper is to find the widest generalization of the NP bound by applying it to 
a target r.v. which is neither Normal nor Gamma, and in both Wiener space and Wiener-Poisson space. 
This is worked out in [lO] but the conditions needed to apply the NP bound are quite restrictive (it was 
also carried out only in Wiener space). The conditions we are introducing here are more general, and are 
still wide enough in scope to cover a Z belonging to the Exponential family or the Pearson family. We 
point out that Wiener-Poisson space is more inclusive than Wiener space (which can be identified with a 
subspace of the former). In fact, it includes processes with jumps, and therefore considers Poisson space too 
as a subspace (also by identification). Nevertheless, even if Wiener space is less general, we can apply our 
techniques to a wider class of target r.v. Z than in Wiener-Poisson space (which requires boundedness of 
the second derivative of the solution of the Stein's equation, something not needed in Wiener space) . 



Our main results are the NP bounds on dy, (A, Z) in Wiener space and in Wiener-Poisson space. The 
main result in Wiener space fTheoremll2|) is 

dniX,Z)<m\g, {X)-gx\ 

< k 



E g,{Xf -E g,{Zf 



|E [AG* (A)] - E [ZG., {Z)]\ + |E [g 



The main result in Wiener-Poisson space (Theorem [26]) is 

dn (A, Z) < fc f E |.g, {X)-gx\+^\{ x (DX) 



-DL-^X 



Here, gx ■= E[(Z)A, — £'L^^A)^|A] and gz E.[(DZ,—DL^^Z)^\Z] are random variables defined using 
Malliavin calculus operators, specifically, the Malliavin derivative D and the inverse of the infinitesimal 
generator L of the 0-U semigroup. The definitions parallel each other, but it would be helpful to think of gx 
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as an object belonging exclusively to X, and gz to Z. On the other hand, g*{-) := E[{DZ, —DL^^Z)f^\Z = ■] 
is a function whose support is the support of Z, taking on nonnegative numbers as values. It will depend 
only on the density of Z, and is independent of the structure oi X. As such, it is an object belonging 
solely to Z. In the second term of the first bound above, G* is an antiderivative of g*, provided it exists. 
We can make sense of the NP bounds in the following way: if we want to know how different the laws 
of X and Z are, then we need to know how different (in the sense) gx = K[{DX,—DL~^X)s^\X] and 
g*{X) ~ M[{DZ, —DL^^ Z)^\Z — X] are. In Wiener-Poisson space, we consider in addition how close the 
jump part E[(|a; [DX] |, \ — DL^^X\) i^] is to 0, which makes sense since Z belongs to Wiener space (subspace 
of the Wiener-Poisson space without jumps). 

In our bounds above, fc is a constant that does not depend on X but on Z and the metric we are 
using. For convergence problems, we do not need its specific value since the convergence will follow from the 
convergence of E \g^^ {Xn) — 5x,J to 0. This presupposes we have such a constant k. This constant appears 
as a bound < k) for some function which is related to the solution of the underlying Stein equation. 

In particular, since we have a Stein equation for each Z (the target r.v.), k depends on Z. Finding such a 
bound k is easy when Z is Normal: 5* is constant, and consequently, the Stein equation is simpler. If g^ 
vanishes at a finite endpoint of the support of Z, the challenge now is to find a bound for To the 

best of our knowledge, [10 (Kusuoka and Tudor) presents the first attempt to find such a sup norm bound 
when Z is not Normal. Their result is presented below as Lemma [T] In Theorem |9] we improve their result, 
and this paves the way for the needed bound we stated for the general non-Normal case. 

The paper is organized as follows. In Section[2l we review the operators we need from Malliavin calculus. 
We also define the functions g.^, and G■^, as well as the random variables gx and 52, studying carefully their 
properties (needed in the subsequent sections). Section |3] contains preliminaries on Stein's method. Here 
we find universal bounds on the first and second derivatives of the solution of the general Stein equation. 
Our main result in Wiener space is in Section |4l where we give a tractable upper inequality which is easier 
to compute. We also characterize when the law of X is the same as that of Z. Said result is applied to 
specific cases when g^ is a polynomial and when is a sequence of multiple integrals. As an example, 

we prove the convergence of a bilinear functional of a Gaussian subordinated field to a r.v. by computing 
some moments and showing their convergence to desired values. In Section [5l we extend the main result to 
the more general Wiener-Poisson space. Here, we work out some sufficient conditions for convergence to a 
Normal law and convergence of the fourth moment. 

2 Elements of Malliavin calculus and tools 

For the sake of completeness, we include here a brief survey of the needed Malliavin calculus objects. The 
r.v. (^DF,—DL~^F^^, to be constructed for = Z and for F = X, is a key element that bridges Stein's 
method and Malliavin calculus. D is the Malliavin derivative operator and L is the generator of the Ornstein- 
Uhlenbeck semigroup. 

2.1 Wiener space 

Nualart presents in Chapter 1 of a very good exposition on Malliavin calculus in Wiener space. We 
mention here the elements that we need. Let ^ he a real separable Hilbert space. Assume a probability 
space (r2, J^, P) over which W — {W {h) : h E S)} is an isonormal Gaussian process. By definition, this means 
is a centered Gaussian family such that E [W (hi) W (^12)] = {hi,h2)^. We may also assume that is the 
cr-field generated by W. The white noise case is when = L? (T, B, /x) where (T, B) is a measurable space 
and /X is a cr-finite atomless measure. The Gaussian process W is then characterized by the family of r.v.'s 
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{W{A):AeB [0, T] ,n{A) < 00} where W (A) = W (Ia); we write Wt ^ W (l[o,t]) ■ We can then thmk of 
W as an (il, J^, P) random measure on {T,B). This is caUed the white noise measure based on fi. 

The g*'' Hermite polynomial Hq is given by Hq (x) = (— 1)'^ e"^^/^^ (^e^^'^^^ for q > I and Hq (x) = 

1. The q^^ Wiener chaos Hq is defined as the subspace of (il) — {ft,T,P) generated by the r.v.'s 
{Hq {W {h)) : /i e ^3, \\h\\^ — l}. In the white noise case ^ = L^^ ([0, 1]), each Wiener chaos consists of iterated 
multiple (Wiener) integrals 

Ig{f)-=n^- I [ ■■■ r ' f{tl,t2,...,tq)dWt^---dWt,dWt, 

Jo Jo Jo 

with respect to W, where / € io®' is a symmetric nonrandom kernel. When / is nonsymmetric, we let / 
denote its symmetrization, and Iq{f) — Iq{f)- 

All elements of Hi are Gaussian and all elements of Ho are deterministic. It is well-known that (il) 
can be decomposed into an infinite orthogonal sum of the Wiener chaoses, i.e. ($7) — (S'^oHq. In the 
white noise case, any F £ (57) admits a Wiener chaos decomposition of multiple integrals 

CO 

where each symmetric fq G Sj®'^ — (T'^) is uniquely determined by F. Note that Iq (/q) — fo —E [F] and 

E[Iqifq)]=OioVq>l. 

Consider an orthonormal system {cfc : fc > 1} in .$3. For / € jo®^ and g e io®"^, the contraction of order 
r < min {p, q} is the element f ®r 9 & ^^®{p+q~2r) defined by 

00 

f'^rg= E (/, (g) • • • «) CiJ^g, (5,eji (g) • • • (g) e^J^^, . 

■il ,...,v 

Even if / and g are symmetric, / (Sir 9 may be nonsymmetric so we denote its symmetrization by f<§r9- In 
the white noise case Sj = L'^ (T), the contraction is given by integrating out r variables. Thus, if / G {T^) 
and qeLl (T«), we have / (gr 5 e (yp+g-2r^ 

if glr 5) (^li ■ • ■ : tp+q-2r) = / f {h,..., tp^r, Si,...,Sr)9 (^p+1, • • • , tp+q-r, Si, . . . , Sr) dfl (si) ■ ■ ■ dfl {Sr) ■ 

Jt'- 

The product of two multiple integrals is 

I, if) Ip {9) = E ^' \r) ®^ 3) ■ (2) 



The Malliavin derivative of a random variable F £ (17) is an ij-valued random variable denoted by 
DF. In the white noise case Sj ^ Ll (T), if F = h (/) = /j, / (t) dWt, then D maps F to an (r)-valued 
element: DrF — f (?') for r £ T. In general, if F e (17) admits the decomposition ((!]), then 



00 



DrF=\qIq^iifqir,-)). (3) 



q = l 



We denote by D^'^ the domain of D in (17). F with the above decomposition is in D^^^ if and only if 



E 



\DF\\%m 



— J2^i 1 ■ 9- ll/gllis (j^g) < 00. D satisfies the chain rule formula: D (/ (F)) = /' (F) 
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when F G D'^'^ and / is differentiable with bounded derivative. One may relax this to almost everywhere 
differentiability of / as long as F has an absolutely continuous law. 

D has an adjoint, the divergence operator S, so that if G DomS C then 6 {F) G {fl) and 

E [i5 (F) G] =E [{F, DG)^] for any G G D^^^. In the white noise case, S is called the Skorohod integral: for 
F G Dom(5 C i^xP ^ ^) ^itl^ chaos representation F (t) — J^'^o^i ifi ')) where each fq G L^g,(<j+i) is 

2 

< oo, i.e. F G Dom(5. 

One other operator we need, L, acts on as in ([1]) in this way: LF — — 9^9 (/"j)- domain 

consists of F for which 9^ ' 9- II/9IIl2(T9) ^ ^ ^-'^^'^ happens to be the infinitesimal generator 

of the Ornstein-Uhlcnbeck semigroup Tt, defined by TfF — X]^o (/"j)- ^^'^ important relation is 

SDF = —LF. More than L, we need its pseudo-inverse defined by L^^F = — q^q ifq)- It easily 

follows that L-^LF ^ F - E[F]. 

2.2 Wiener-Poisson space 

Assume a complete probability space (il, F, P) over which C = {£t}oQ is a Levy process. By definition, this 
means C has stationary and independent increments, is continuous in probability, and Cq = 0. Suppose C 
is cadlag, centered, and E < oo. We may also assume F is generated by C. Let C have Levy triplet 
(0,cr^,i^) and thus, Levy-Ito decomposition d — o-Wt + J J^^ ^^^^^xdN {s,x) where W — {Wt}(>Q is a 

standard Brownian motion, N is the compensated jump measure (defined in terms of v) and Rq = M — {0}. 
See [U and [23] for more about Levy processes. 

Consider now the measure fi on B (M+ x M) where M+ ~ {t : t > 0} and 

dfi {t, x) ~ a^dtSo (x) + x^dtdv [x) (1 — (5o (x)) . 

Analogous to a Gaussian process W being extended to a random measure (which we also denoted by W) in 
Wiener space, C can be extended to a random measure M (see 9 ) on (M+ x M,f?(M+ x R)). This is used 
to construct (in an analogous way to the Ito integral construction) an integral on step functions, and then 
by linearity and continuity, extended to L^g,, = ,B (R+ x R)'' , ^«"?) . We also denote it by Iq. 

As in Wiener space, 

1. /,(/)=/,(/); 

2. /g is linear; 

3. E [Iq (/) (.9)] - 

Thus, when F = /,(/), E[F2] = Wqiff] = Q^WM^,.- 

Contractions are defined slightly differently. Suppose / G i^«<, and g G ^^®p. Let r < miii {q,p} and 
s < min{g,p} — r. The contraction / (g)^ 5 G L^^(,+p_2r_s) is defined by integrating out r variables and 
sharing s of the remaining variables: 

U ®lg){z,u,v) = Wx^ if {■,z,u) ,g{-,z,v))^2 



symmetric in the last q variables, 6 (F) = J2q=Q ^9+1 [fq] if J2q=Q + fq 
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where z e (IR+ x Mq)^, Zi — {ti,Xi), u £ (R+ x RqY * and v G (M+ x Mq)^ Its symmctrization is 

f®1.g. We need the following product formula later (see [TT] for the proof): 

phqpf\q-r / \ / \ / \ / \ 

I. if) I. (9) = E E/!^! (') 7) : Vp-..-. (/ K 9) ■ (4) 

We may think of this as a more general version of the product formula ^ where we only consider s = 
since there are no jump components to be shared (which appear in the definition of / (g)^ g). 

We have briefly narrated a setup parallel to what was done in Wiener space. See [24j for a more detailed 
exposition. This time though, we have only considered S) ~ (M+ x M, S (M+ xM.) , ^) as underlying 
Hilbert space, with inner product {f,g)^^ = Jjg_i_ f (z) g (z) dfi (z) . There is as yet no Malliavin calculus 
theory developed for a more general abstract Hilbert space. While we don't have a chaos decomposition 
via orthogonal polynomials (like Hermite polynomials in Wiener space; see [7]), we still have a comparable 
decomposition proved by Ito (Theorem 2, [9 ): for F £ (57, J^, P), 

oo 

^ = E ^9 (/g) where fg £ i^^, . (5) 

q=0 

With this decomposition, we can define the Malliavin derivative operator and Skorohod integral operator. 
Define DomI? as the set oi F £ (ft) for which J^^i 11^- < ^ and 

oo 

i?,F = ^qVi(/,(z,.)). 

9=1 

It is instructive to consider the derivatives Dt^ and Dz where z = (i, a;) has x ^ Q. This will enable 
us to better understand the similarities, and where they end, between the Malliavin calculus of Wiener 
space and that of Wiener-Poisson space. See [24] and [25] for more details on the following discussion. 
We consider two spaces on which we can embed DomZ). For F £ L'^{VL), we say F £ DovaD'^ iff 
Er=i99!/K+IIA((^.0)>-)lli. <oo andFeDomi^-^ '^^ H^^M U .wj U ^ < 

OO. In fact, DomZ? = DomD" n DomZ?''. Since W and N are independent, we can think of f2 as a cross 
product of the form flw x flj where ilw — C (M+) and fij consists of the sequences {{ti.Xi) , (^2, X2) , ■ ■ ■) £ 
(R+ X Mq) (with a few other technical conditions). 

• The derivative Dtfl can be interpreted as the derivative with respect to the Brownian motion part. In 
fact, if = 0, then Dt.oF = F where is the classical Malliavin derivative (defined in Wiener 
space); the i comes from the fact that we are differentiating with respect to aWt and not just Wt- 
From the isometry (ft) ~ (^w'l (^j)) , consider F £ L'^ (fl) as an element of (^w'-: (^j)) ■ 
A smooth F then has form F — GiHi where each Gi is a smooth Brownian random variable 

and H, £ {flj). We can then define by D^F = J27=i {D^G,) iJ„ where D^G, is the 
classical Malliavin derivative. It can be shown that this definition can be extended to a subspace 
DomD*^ C DomZ?°, so that for F £ DomZ?'^, as expected, 

Dt^oF = -Df'F. 
a 

For functional of the form F ^ f {G,H) £ (f^) having G £ DomD^, H £ L'^ {^j), and such that 
/ is continuously differentiable with bounded partial derivatives in the first variable, we have a chain 
rule result: F £ DomD" and Dt.oF = jrfj {G,H) D^G. We may loosen the restriction on / to a.e. 
differentiability if G is absolutely continuous. 



6 



The derivative Dz, z ~ {t, x) with a; 7^ 0, is a difference operator: for F E DomD^^ 

^ ^ _ F jujt,,) - F juj) 



where, if "^zF is the right-hand expression, then E /r+xEq dij{z) < 00. The idea is to 

introduce a jump of size x at time t which is captured by the realization ujt,x- For oj = (w^,a;''), 

we define cjt.x by simply adding the time-jump pair {t,x) to uj'-' . For F — f {G,H) g (ft) with 
G E (^j)i H E DomD'^ and / continuous, we have this chain rule result: 

jj^p ^ nG,H{ut,x))-f{G,H{u)) ^ f{G,xDzH + H{uj))-f{G,H{oj)) 

X X 

If / is differentiable, then by the mean value theorem, for some random 6^ E (0, 1), 

df 

D,F = (G, e^xDzH + H {u)) D.H. 

The following unified chain rule will be very useful (see Proposition 2 in [3D]): If F e Dom Z)^ nDom Z?-', 
DF E L'j^, f E C''^^ has bounded first derivative (or /' may be unbounded if F is absolutely continuous with 
respect to Lebesgue measure) and f^^~^^ is Lipschitz, then for z E {t, x) E K+ x M, 



Dzf{F)^y^^-^x^-\DzFr' I ^ ^ 



- XU] 



(fc-1)! 



In the case where /('^ is differentiable everywhere, the chain rule is 



n=l 



fc! 



(6) 



(7) 



for some function 0^ E (0, 1) for all z ~ {t, x) E 



We now define the adjoint of D (see [21] again). Suppose F E L'^ {M.+ x M. x ft, B {R+ x R) x F, n x P) 
with F (z) = X](^o ^9 (/<? (-^j ')) where each fq E i^g,;^^!) is symmetric in the last q variables. In this case, the 

< 00, i.e. F E T)omd 



Skorohod integral of F is 5 (F) = J2T=a h+i ( fq ) where Y^n^a il + 1)' fq 



(by definition). Furthermore, E [6 (F) G]^E {F, DG) ^2 



for any G E DoniD. 



Finally, we define as before L = —SD: for F as in (O, LF = — X^^i l^q ifq)- The pseudo-inverse is 
defined by L^^F = - J2'^^i \lq ifq)- We have again L-^LF = F - E [F]. 



Remark 1 Write Zq — (zi, . . . , Zq), with Zi — {ti, Xi) for all i. Define 

W — <F = ^ Iqifq) e DomF)" : fq E L'^^e>q, and for every q, fq{zq) — if Xi ^ for some i > . 
I q=0 J 

Notice from the previous discussion that if fq{zq) — because Xi 7^ 0, then Iq{fq) coincides with an iterated 
multiple (Wiener) integral. Therefore, Wiener space can he seen as a subspace of Wiener-Poisson space 
(similarly for Foisson space as a subspace). Moreover, W coincides with the subspace D^'^ (through embed- 
ding). The relevance of these facts is that if we have a r.v. F E W , then the chain rule formula and the 
Malliavin calculus operators are exactly (up to a constant) the same as those in Wiener space (as explained 
earlier in this subsection). Furthermore, the results (from other papers) in Wiener space can be replicated in 
W and so the conclusions will hold in Wiener-Poisson space, but within W . From now on, D^'^ will mean 
the subspace D^'^ in Wiener space or the respective embedding W in Wiener-Poisson space. 
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2.3 The random variables gx, Qz and the functions G* 

From this point on, f) will be taken as I? (K+ x R, (M+ x M) , /i) if we are in Wiener-Poisson space. Now 
suppose F has mean 0. We have the following integration by parts formulas. 

• If F e DomD^ n DomZ?'' and f E with Lipschitz first derivative assumed to be bounded (or /' 



may be unbounded if F has a density), 
E [Ff (F)] = E \{-DL-'F, DF)^ f {F) 



+ E 









DF 



i5- 



(8) 



• If F e DomZ)'^ n DomZ)'^ and / is twice differentiable with bounded first derivative (or /' may be 
unbounded if F has a density), 

f"{F + 0.xDF) 



E[F/(F)]=E {-DL-'F,DF) f{F) 



■E 



-DL-'F, 



-x{DFy 



(9) 



If F € D"'^'^ (see Remark [T]) and / is differentiable with bounded derivative (or / is at least a.e. 
differentiable if F has a density). 



E[F/(F)]=E {-DL~^F,DF) f{F) 



(10) 



These formulas provide the link to using Malliavin calculus techniques in solving problems related to 
Stein's method. Since F = LL^^F = —SDL^^F, we have 



E [Ff (F)] = E [SDL-^F ■ f (F)] = E {-DL-^F, Df (F)) 

A direct application of the chain rule for Wiener-Poisson space, choosing k — 2 in ^ and ([7]) , yields 
(HI) respectively, and an application of the respective chain rule in Wiener space yields pO)) . 



and 



Recall that in the Wiener-Poisson case (by Remark [T]), if F e D^'^ then the chain rule formula ([7]) will 
reduce to the corresponding one in Wiener space. Note that in this paper, we are assuming the target r.v. 
Z is in D^'^. Consequently, the points we will make about Z will be valid whether we are working in Wiener 
space or Wiener-Poisson space. The discussion of this subection then applies to both spaces. 

Assumption A Z E D^'^ has mean and support {l,u) with — oo < I < Q < u < oo. The density p* 
of Z is known, and it is continuous in its support. X is either in D^'^ (Wiener space case) or in 
DomZ?'^ n DomZ)"' (Wiener-Poisson space case), and it also has mean 0. 

Caution: Notice that in the previous subsection we used a; e M to denote the jump component of z e x K 
in our state space. On the other hand, we are using Z to denote the target r.v. and X the r.v. with unknown 
distribution. A confusion may arise in the usage of x and X, or z and Z. However, we will stick with current 
notation for consistency with existing literature. In this regard, we urge the reader to keep in mind that x 
represents the size of the jump while AT is a random variable not (directly) related to x. On the other hand, 
z is a jump (time of the jump, size of the jump) while Z is the target r.v. which has no jumps. 

Remark 2 In some results, we will consider instead of X a sequence {A„} of random variables. In this 
case, we have the same assumptions (and corresponding functionals, defined below) for each A„. Note that 
for Z G D^'^, the support necessarily has to be an interval (see Theorem 3.1 U8f . Proposition 2.1.7 JWl), 
a consequence that carries over to Wiener-Poisson space. The continuity assumption of the density is 
not strong at all, since general processes like solutions of stochastic differential equations driven by Brow- 
nian motion or (under mild conditions) fractional Brownian motion (for example see 121) have continuous 
densities. 
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Define the random variable gp — ¥, {DF, —DL ^F)^^ F for any Malliavin differentiable r.v. F (we 
win work with F = Z and F ^ X later) . Nourdin and Peccati proved that > almost surely (Proposition 
3.9, [14]). Closely related is the function 

5*(z):=E {DZ,-DL-^Z)^ Z^z . (11) 

Trivially, gz — g* (Z). Nourdin and Viens (Theorem 3.1 181) proved that Z has a density if and only if 
{Z) > a.s. Therefore, (Z) > a.s. (Assumption A) and 5* (z) > for a.e. z £ {l,u). Equation (ITU)) 
implies E [Zf (Z)] = E [gzf (Z)]. In the same manner, in Wiener space {X e D^^^) we have E [Xf (X)] = 
E[gxf' (X)], while in Wiener-Poisson space {X e DomD^ n DomD''), this changes to E[Xf{X)] = 
E [gxf (X)] + E[{-DL-'^X, ^" ^^^l'"''^^'' x{DXf) h]- One needs to be careful not to write E [Xf (X)] = 
E [g^ {X) /' {X)] or 5* (X) ~ gx a.s., both false, since is derived from the law of Z (we would need the 
corresponding g^, of X to make it true). 

Nourdin and Viens proved (same Theorem 3.1 [T8^) that 



VP* iv) dy / yp* (y) dy 

g*{z)^^ = fora.e. ze (Z,u). (12) 

p* [z) (z) 

In their proof, they pointed out that <^{z) := yp^ {y) dy — ~ yp^ (y) dy > for all z £ {I, u). Since 
is (necessarily) bounded (Assumption A), {p{z)/ p^{z) is strictly positive (inside the support). From this 
point on, we will take 5* to be either (jlll) or the version (jl2p . whichever suits our purposes. Furthermore, 
we can assume that g* (z) > for every (and not just for almost every) z £ {l,u). Notice that using this 
definition of 5* we can conclude that ((7*(z)p*(z)) = <^'{z) = —zp^,{z). 



Given the density of Z, we can compute using (|12p . Some examples of known distributions with 
their (7* are given in Table [TJ Recall that g*{z) = outside the support. 

Conversely, one can retrieve the density given g^, using the following noteworthy density formula 
Nourdin and Viens [T^ proved: 

E\Z\ ( r y \ 

P* = — TT f^xp - / — i—rdy . (13) 



2g* {z) \ Jo g^ (y) 
Proposition 3 necessarily satisfies the following: 

/ — T~^dy = -00 / —^dy = 00. (14) 
Ji 9* (y) Jo 9* (y) 

Proof. With (p{z) defined as before, Nourdin and Viens (Theorem 3.1 ^8j) showed that 

y -dy = ln^^'^ 



9*iy) ^(z)' 
Since iy9(z)— >-Oasz— >m and as z — !• the result follows. ■ 

Remark 4 The necessary conditions in Proposition^^ are actually not new. Stein (Lemma VI. 3 !26J ) has 
pointed out that these are necessary for !il2\} to hold. 



Suppose g^ {x) — a {x — lY for some constant a > and the support of Z is (/, 00). Then -^-j^dx — 

00 if and only if p < 2, and J^ ^ ""■'^ ^''^^V */ ^ ^P- Similarly, if g* (x) — a {u ~ xY over 

the support (— oo,u), 1 < (7 < 2 necessarily. Also, if g* {x) = O {x^) and the support is (—00,00) , then 
P !i 2. If g* (x) — a{u — x^ {x — lY over the support {I, u), then p > 1 and q>l necessarily. 
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• Not every satisfying Jj^P will belong to a random variable in D^'^. For instance, suppose Z is 
Inverse Gamma (see Table{l^ having support (l,oo). For this random variable, g^ (z) — a{z ^ I) for 
some a > 0. Z can be shown to have finite variance if and only if a < 1. Thus, when a > 1, is 
satisfied but Z ^ L'^ (57) so Z ^ D^'^. This means that the coverage of our method is not as extensive 
as we would like it to be. Hence, it is important that we check if the target r.v. belongs to B""^'^. 



Z , with support and parameters 



9* (z) 



Normal 

(— oo, oo): cr > 



Gamma 

{I, oo): I = —rs, r > 0, s > 



1 , ,\r-l I z — I 
(z — t) exp 



sT (r) 



s{z-l) 



X 

{I, oo): I = -V, d.f. V > 



2''/2r(i>/2) 



(z — Z) ^ exp 



z-l 



2{z-l) 



Exponential 

{l,oo): I = -i A > 
A 



Aexp (—A (z — I)) 



Beta 

{l,u): I = - 



r + s 



u = l + l,r, s>0 



13 (r, ,s) 



{z-iy-'{i + i-zY 



{z-l){l + l-z) 



Pearson Type IV 



C(l + {z~tf^ '^gsarctan(2:-t) 



^(l + (-t)^) 



2(r-l) 



Student's T 

(— oo, oo): d.f. V > 2 



(!) 



^ 1 + ^ 

w — 1 V V 



Inverse Gamma 
{I, oo): I = — 



r-1' 



r > 3, s > 



S' r ( S 

(z — I) exp — - 



r(r-i) 



' {z-lf 



r-2 



Uniform 

{l,u): u = -I > 



1 

2u 



Pareto 

(;,oo): c> 2, / < 



(z - cl) 



c+1 



c - 1 



(z -l){z- cl) 



Laplace 

(— oo, oo): c > 



■exp(-c|z|) 



(l + c|z|) 



Lognormal 

(;,oo): l^ - exp{S+ \a'^) 



r,=—T7 exp ( [p (z) + af 

In (z - /) - (5 
where p \z) ~ 



ae^^ exp 



X / e I'^ds 

lp{z)-(y 



Table 1: Common random variables Z with their p^^ and g■^, 



Parallel to g^ of Z , we may define a corresponding object for X but we will have no use for it. In fact, we 
typically won't have access to the density of X; if it was known otherwise, one may characterize X without 
having to approximate it by Z . We will only assume properties (for X) amenable to the Malliavin calculus 
that would allow us to define gx- 
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Let G* (z) = /;^.9* (y) dy be the indefinite integral of 5* (assuming € L^{l,u)). Consider the Wiener 
space case {X G D^^^) and suppose lls*!!^^ < 00 or X has a density. If we take / = G* in (ITUl) . then 

E[gx5* W] =E[a(X)X]. (15) 

Assumption A' Along with Assumption A, either < 00 or A has a density. 

Proposition 5 (Moments formula) 

• Wiener space: E = rE [F^^^gp], provided the expectations exist. 

1. If g^{Z) ^ gz is a polynomial in Z , i.e. 5*(z) = OfcZ*^, t/ien E = X^fcLo '^'^fe-'^ [^''^'''~'"] ■ 



2. If X = Iq {g), then E [X''+^] 



A'-i IIDAII?. 



Wiener- Poisson space: E = rE -'-(7^?] 



r(r — 1) 



E [{-DL-^F,x{F + e.xDFy''^{DF)'^)s,] 



1. If X ^ Iq [g), then E [A''+i] = i^E 



A"-! IIZJAI 



'•('•-i^E 



x{DXf ,{X + e.xDX) 



Proof. Simply let / (y) = in ([9]) and (jlOp. taking into account the fact that —DL ^F — -DF when 
F^IM- ■ 



3 Stein's method and the Stein equation 

Stein's method is a set of procedures that is often used to measure distances between random variables such 
as X and Z . More precisely, we're measuring the distance between the laws of A and Z . These distances 
take the form 

d-H (A, Z) = sup |E [h (A)] - E [/i (Z)] | (16) 
hen 

where 'H is a suitable family of functions. If we take T-Lw = {h ■ ^ 1} where ||-||^ is the Lipschitz 
seminorm, then dw = d-Uw is called Wasserstein distance. The bounded Wasserstein (Fortet-Mourier) 
distance corresponds to V-fm = {h : ||^||i + ||^|loo — Clearly, dpM < dw- dpM is important because 
it metrizes convergence in distribution: dpM {Xn,Z) ^ if and only if A„ ^ Z. dw on the other hand 
induces a topology stronger than that of convergence in distribution. 

Nourdin and Peccati [1^ mentioned other useful metrics. We have the Total Variation distance when 
'Htv — {1b '■ B is Borel} and the Kolmogorov distance when Hk — {l(-oo,z] '■ z G M}. The latter for 
example is suited for the analysis of probability tails. However, in this paper, we will only consider dw and 
dpM as we try to find bounds for d-H (A, Z) by exploiting properties of Lipschitz functions h G H. 

A Stein equation is at the root of Stein's method. Given Z and a test function h, the Stein equation is 
the differential equation 

g, ix) f (x) - xf {x) ^h{x)-¥.[h {Z)] (17) 

having solution / — fh- If the law of A is "close" to the law of Z, then we expect E [h (A)] — E [h [Z)] to 
be close to 0, for h belonging to a large class of functions. Consequently, E [g* (A) /' (A) — A/ (A)] would 
have to be close to 0. In fact, subject to certain technical conditions, the left-hand side of equation ([T7|) 
provides a characterization of the law of Z: E [g* (A) /' (A) — Xf (A)] = if and only if A Z (in the 
equation, information about the law of Z is coded in g.^,). The following proposition states this result in its 
precise form. For a quick proof, see Proposition 6.4 in }14 . The first statement is Lemma 1 in 26 by Stein. 
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Lemma 6 (Stein's Lemma) 

1. If f is continuous, piecewise continuously differentiable, and E [g* (Z) |/' (Z)|] < oo, then 

E[.g* {Z)f'{Z)-Zf(Z)] = Q. 

2. If for every differentiable f , x ^ jg* {x) /' {x) \ + \xf {x)\ is bounded and 

E[g, iX)f' {X)-Xf{X)]=0, 

then X Z . 



(18) 



(19) 



Let v. = Hfm or % = 'Hw Using (fT7|) on pBj) . we have 



{X, Z) < sup |E [.g* (X) f (X) - Xf (X)] \ 



(20) 



where the sup is taken over the family 3^-^ of all Stein equation solutions / corresponding to h Cz H. Here the 
integration by parts formulas ([8]) and pO|) allow us to rewrite the term E[X/(X)] in terms of the derivatives 
of / and the r.v. gx, as we pointed out before. For instance, in Wiener space, 



dn (X, Z) < sup |E [g, (X) f (X) - gxf (X)] \ = sup |E [/' {X) [g, {X) - gx) 



(21) 



Thus, to ensure that the distance between X and Z is small, 5* {X) should be close to gx- We also need 
to have a good control of f (X). One way of addressing this, taking note of Corollary 6.5 in [T3], is by 
assuming a universal bound for E /' (X)'^ for all / e J'-u since 



dn {X, Z) < 



I sup E 

fe'Ju 



f (xy 



E 



{g, (X) - gxY 



(22) 



The first factor is intractable since it requires us to consider conditions on X in relation to all members 
/ of the family 3^-^. If however we have a uniform bound for /', then we can avoid imposing an additional 
restriction on X. In this case, we only need worry about how close {X) is to gx in (ft). In fact, such 
a bound allows us to just consider how close (X) is to gx in (fl). It is then interesting to see how 
information about the law of Z is contained in its Malliavin derivative. Questions of how close the law of X is 



X 



is to E 



{DX,-DL-^X)^ 



X 



Notice 



to that of Z is passed on to how close E {^DZ,—DL ^Z)^ 
though that this discussion needs to be modified slightly in Wiener-Poisson space, since the integration by 
parts formula ([8]) involves also the second derivative. Thus, we need to control (in a uniform way) both the 
first and second derivatives of the solution of the Stein equation. Due to this extra requirement, as will be 
seen later, we will not be able to apply our tools to as wide a scope of target r.v. Z, as we would be able to 
do in Wiener space. 



3.1 Bound for /' 

The Normal case in Wiener space: 

If Z is standard Normal (g*{z) — 1), the Stein equation is f {x) — xf (x) — h{x) — E [/i (Z)] and it has 
solution f{x) = e'''' /"^ jZ^[h{y) - ¥.[h{Z)]]e-y'' /'^dy. Stein proved (Lemma II. 3 in JMI) that \\f'\\^ < 
2\\h-E[h iZ)]\\^. In faci, \\f'\\^ < min {2\\h-E[h (Z)]||^ , 4 \\h'\\^} (see Lemma 2.3 JEI)- For h G Hfm, 



l/'ll^ < 4. It follows from fHp that dpM {X, Z) < kE [|1 - gx\] < E (1 - gx) with /c = 4. Similar 
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estimates for h G Hw ieo-d to a bound for dw of the same form hut with k = I (Lemma 4-2 J^, Lemma L2 
114V - How close the law of X is to the standard Normal law depends on how close gx is to gz = 1 (in the 
L^ sense). 



In the general case, the Stein equation ([TTl) has solution 



1 



[h (y) - ruh] p* (y) dy 



[h (y) - mh] p* (y) dy (23) 



5* (x) {x) J I 



.g* {x) (x) 



X 



for X € (Z, w), where ruh '■= E [h (Z)]. 

The proof of the bound for /' when Z is Normal can be adapted to find a constant bound for g^, f in the 
non-Normal case. If g^, is uniformly bounded below by a positive number, we easily get a uniform bound 
for /'. Unfortunately, this is not always the case. In Table[T]we can see several examples of target r.v.'s for 
which (7* can get arbitrarily close to in its support (for example, when Z is Gamma and (z) — s {z — l)^). 
Kusuoka and Tudor in [TU] (Proposition 3) proved the following proposition to address this issue. We state 
it in the following form using notation and assumptions we have set. 

Lemma 7 Suppose we have the following conditions on <?* . 

L If u < oo, then lim^_,^^ g^ (x) / {u — x) > 0. 

2. If I > — oo, then lim ^^; g^ [x) / {x ~ I) > 0. 

3. If u ~ oo, then lim ^,_,^, g^ (x) > 0. 

4- If I = — oo, then lim .^_,; g^ (x) > 0. 

Then the solution f of the Stein equation JiTp , for a given test function h with \\h\\^ < oo and \\h'\\^ < oo, 
has derivative bounded as follows: 



where the constant k depends on Z alone, and not on h. 

Unfortunately, conditions 1 and 2 are too restrictive. Consider for instance a r.v. Z with support 
(Z, oo) and g^, (x) = a {x) (x — l)^ , where a{x) is uniformly bounded below by some ao > 0. From Remark 
El 1 < p < 2 necessarily. Among all of this form. Lemma [7] is thus only able to assure the needed 
boundedness of /' when p = 1. For instance, when Z is Inverse Gamma or Lognormal, condition 2 fails 
(see the corresponding g^ in Table [1]) . This stresses the need for less restrictive conditions on g^ that would 
allow us to include these cases and much more. The first requirement in order to achieve this is a good 
representation of the derivative /'. 

Proposition 8 For x G {l,u), the derivative f of the solution of the Stein equation |J7[ ) is 



||/'IL<Mll/ilL + ll^'IL) 



(24) 



/' (^) 



5^ (a;) (x) 



1 




[1 - $ (s)] $ (t) [h' (t) - h' (s)] dtds 



where ^{x) = p^,{t)dt is the cumulative distribution function of Z . 
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Proof. First, 



h (x) — nih = [h{x) — h (s)] (s) ds + [h (x) — h (s)] (s) ds 



h' (t) dt 



h' (t) dt 



p* [s) ds - 

pX r ft 1 nU r nU 

= / p* (s) ds h! {t) dt - I IP* (s) ds 

Jl iJl J Jx iJt 

= ^ (t) h' (t) dt - J [l-^{t)]h' {t)dt 



p* (s) ds 
h' {t) dt 



and so 



g* {x) p* (a;) f{x)=j^ [h {y) - m^] {y) dy 



$ {t) h' (t) dt 



P* (y) dy ■ 



{t)dt 



P* (y) dy 



$ {t) h' (t) dt 



/X r pt 1 pu r px 

/ p.{y)dy [l-'^{t)]h'{t)dt- / p,{y)dy 
, Jl J Jx Ui 



P* (y) dy 



[1 - $ (t)] h' (t) dt 



nX pX pU 

= / [^{x)-^(t)]^{t)h' {t)dt- / <^{t)[l-^{t)]h' {t)dt- / ^{x)[l-^{t)]h' {t)dt. 

Jl Jl Jx 



Cancelling some terms and solving for / 
Observe that if x < 0, 



g» (x) p* (x) J I 



$ (t) h' {t) dt - 



$(x) 



g» (x) p, (x) 



[1 - $ (i)] /i' (t) dt. 



(25) 



while if a: > 0, 



pU 

tp4t)dt + / tp,(t)di < x$(x) + 5(*(x)p*(x) 



= E[Z]= tp^{t)dt+ tp*(t)(ii > -g*(x)p*(x) +x[l - $(x)]. 



Therefore, < — x$ (x) < (/«(x)p*(x) ^ as x ^ Z and < x [1 — 3> (x)] < g*(x)p*(x) — ^ as x — ^ w. When 
we then integrate by parts, 

/•X X i^x 

I ^{t)dt = t^{t) - / ip* (t) = x$ (x) + 5, (x) p, (x) (26) 
Jl I Jl 

1 [l-<i{t)]dt = t[l-^{t)] + / tp,{t)dt = -x[l-^{x)]+g^{x)p^{x). (27) 

J X X Jx 



Finally, 

5« (x) /' (x) = xf (x) + /i (x) - m/i 
x[l-$(x)] 



$ (t) /i' {t) dt ■ 



x$ (x) 



+ 1 



r[i-^it)]h' 

J X 



(t) dt 



g* (x) p* (x) J Jl Vfl* (a;) P* (a;) 

= (J [l-^{s)]ds <^{t)h'{t)dt- ^{t)dt [l-^{s)]h'{s)ds 

g» (x) p, (x) J I g» (x) p, (x) J, 
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which leads to the given form of /'. ■ 



The bound (|24p is not directly suited for dw where we don't have a prescribed bound on ||/i||oo- A 
workaround, as pointed out in [10 , is that for each h € "Hwi we pass on the analysis to a sequence {/in} 
converging to h uniformly in every compact set, where {/i„} C {/i G Cq : ll^i'lloo — l}- However, with the 
help of the previous lemma, we can overcome this complication by giving a bound for /' in terms of only 
||/i'||oo- Recall that if h is Lipschitz, it is a.e. difFerentiable and ||/i'||oo — II^'-IIl- Thus, the upper bound 
obtained here is immediately well suited for all / e 'Jfm and for all / € 7w- 

Theorem 9 // applicable, assume conditions 3 and 4 from Lemma Suppose there exists a positive 
function g G C^{l,u) such that 

1. < lim^_j.^ (x) /g (x) < liuix^u g* (x) /g (x) < oo and g'{u^) := lini2._5.„- g' {x) € R exists^ 

2. Q < lim^_^j g^f [x) jg (x) < lima;_j.i g^, (x) jg (x) < oo and '■— limj._>./+ g'{x) € R exists. 

Then the solution f of the Stein equation ^7^ , for a given test function h with \\h'\\^ < oo, has derivative 
bounded as follows: 

ll/'lloo<fc||/^'lloo (28) 

where the constant k depends on Z alone, and not on h. 
Proof. First note that from Proposition [U 

Fix I' and u' s.t. I < V < Q < u' < u. Since g^,{x)p^,{x) is continuous and strictly positive on 
it attains its minimum ni := inf [;/ „/] 5* (a;) > on this compact set. Also by continuity of the 
density M sup[,,,„,] (x) < 00, and .g, (x) = ^^^^ > f > on [I' ,u'], so gl{x)p,{x) > 2^. 
By the continuity and positivity of Ii{x) :— /J [1 — $ (s)] ds and l2{x) := ^ (t) dt we conclude that 
K sup[,,„,] ihix) V hix)) < 00. By m, l/'(^)l < ^ II^'IL on [l',u']_ 

Since I' and u' were arbitrarily chosen, we only need to prove now that lim^:-!.; |/' (a;)| < fci||/i'||oo and 
Imix^u I/' (a;)| < fc2||/i'||oo for some finite constants fci and k2. Due to the symmetry of the arguments it 
suffices to prove just one of these limits. Suppose I' was chosen small enough so that g & {I, I'), and for 
some constants < c < C < 00, cg^ (x) < g {x) < Cg^ (x) on {I, I'). 

• Case 1: I > —00. 

We show that the limit of the right-hand side of ([29]) is finite as a; — >■ I. Note that in this case, 
[1 - $ (s)] ds = (x) (x) ~ x[l - ^ (x)] \l\. By L'Hopital's rule, 

Ih^ I/' (x)| < 2 \\h'\\^ \l\ li^ 5^/^ ^(^)dt ^ 2 ^ $ (x) 



*i °° x^i g (x) (x) (x) °° x-^i —xg {x) [x) + g' {x) g* {x) p* {x) 

x^i [—cx + g [x)\ g^ (x) (x) g [i^j — cl x^i — .xp* [x] g ((^j — ct 



CO 



Since g{l^) '■= lim2_>.;+ 5* (2;) = and 5 > 0, we may assume V is small enough so > on {1,1'). 
Consequently, g' ^ cl < 0. 



stands for the extended real numbers, i.e. R - 
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• Case 2: I — —oc. 

Since lim ^ , ^ g^, (x) > 0, we may suppose I' is small enough so that for some constant mo > 0, 
g* (x) > mo over (— oo,/'). 

(g* (x) (x) - X [1 - $ (x)]) r $ (t) dt 

°° \x^-oo mo x^-oo gi[x)p^[x) 

\x\ r <^>{t)dt 

x^-oo gi{x)p^[x) 

There are two subcases to consider depending on the behavior of 5(2:) as x ^> —00. From the continuity 
of g and the existence of 5' (/"*"), L := lima; -i.-oo 5(2^) necessarily exists. If L < 00, then lim3;_j._oo = 
0. If L = 00, then by L'Hopital's rule, lim^^_oo ^l§^ = — ^^'<^x^-oo g' (x) — —g'{l^). In either case, 



lim^r^-oo ^ exists. 

— Subcase 1: lim^_s._oo — 00 



Note that by (gHl), <^{t)dt = x^{x) + g^,{x)p^,{x) < g^{x)p^{x) so 

1^1 Too ^it)dt ^ ^ \x\g,{x)p.{x) ^ ^\x\ 



Therefore 



gl{x)p*{x) g{x)g^{x)p^{x) g{x)' 



lim \nx)\<2\\h'\\^C hm = < 00. 

:— >■ — 00 x^ — CO g[xj 



Subcase 2: lim^_>._oo 4^ < 00 
Similarly from ([26l) . 



\x\r_^mdt ^ f^g.{t)p.,{t)dt ^ j:^g,{t)p,{t)dt 

gl{x)p*{x) ~ mog*{x)p*{x) ~ mog*{x)p*{x) 

Therefore, 

li^ l/'(.)l<»^ li^ ^l^l;f^<m\^ li^ 



mo x^-ca (7*(x)p*(x) mo a;-j--oo — a;p* (x) 

< hm — — - < 00. 

mo x^-00 c\x\ 

The proof that lim2;_>u |/' (a;)| < fc2||/i'||oo for some < 00 is similar. ■ 

Note that if g* is uniformly bounded below in a neighborhood of ^ > — cxj (or for m < 00) then condition 
2 (1 in the case of u) from Theorem [9] is not required (see discussion before Lemma [7|). In the statement of 
the previous theorem, we can take 5 = 5* if g* is continuously differentiable (at least locally close to the 
endpoints of the support), and in this case the conditions are trivially met. In other words, if we can check 
that 5* € C^{l,u) then bound (|28|) is automatically true (given the existence of g'{u~) and g'(Z+)). These 
new conditions are met by all r.v.'s in the Exponential family, Pearson family, and practically any other r.v. 
whose density is and is strictly positive in its support. If g» is not continuously differentiable, we can 
still get the bound but we are required to approximate 5^ by a continuously differentiable function g near 
the endpoints of the support. For example, consider the Laplace distribution where g^,{x) = ^(1 + c|x|) (see 
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Table [T}. In this case is differentiable everywhere except at 0. Therefore we can choose g{x) — g*{x) for 
all x G {—oo,l') U (u',oo) (with — oo < I' < < u' < oo) and g{x) = <t>{x) on {l',u') where is a smooth 
function such that g is differentiable at V and u'. 

Assumption B We have the folfowing conditions on g^, . 

1. For some positive g € {l,u), 

(a) < lima.^„5* (x) /g (x) < hm^r^u g* (x) /g (x) < oo. 

(b) < lim^^, 5* (x) /g (x) < liuix^i g* (x) /g (x) < co. 

(c) g'{l^) and g'{u~) exist. 

2. If u — oo, then hm ^^,, {x) > 0. 

3. If / = — oo, then hm ,^ q^, (x) > 0. 

3.2 Bound for /" 

For our convergence in distribution resuhs in Wiener-Poisson space, we need a boundedness result for /". 
The existence of /" demands more conditions on g^ such as differentiability, which is understandable since 
we are requiring greater regularity in the solution of the Stein equation. In this setting, the existence of /" 
will also immediately force most conditions of Theorem [9] to be satisfied. If we want to work with dw or 
dpM, we need to consider Lipschitz functions h, and for any such test function, we can only hope for it to 
be differentiable almost everywhere. Consequently, /" must be understood in the almost everywhere sense, 
i.e., /" is a version of the second derivative of / such that wherever the second derivative does not exist, /" 
will have a value of 0. 

Before setting out to find a bound, we point out the unfortunate fact that our results here will not apply 
to as wide a range of target r.v. Z as what happened for the first derivative. More specifically, we won't be 
able to give a finite bound for \ f" [x)\ when I > — oo or u < oo, as we were able to do for |/' (a;)| in Theorem 
[9l We actually have a counterexample to illustrate this: a r.v. Z with support (^,oo) M, such that for 
some Lipschitz and bounded /i, /" {x) does not tend to a finite limit as a; — /. A similar counterexample 
can be constructed for a r.v. Z with support (— oo, u) M, or with support {I, u) ^ M. 

First, we make preliminary computations on /". Differentiating ()17|) gives us the second derivative 

/" {x) - ^^—^f {x) + ^/ (x) + ^h' (x) 
9* [x) g^ (x) 5, (x) 

which, after considering the form of / in equation (|25p and of /' given in Proposition |51 reduces to 



A{x) / ^{t)h' {t)dt + B{x) / [l-^[s)]h' {s)ds + gl{x)p,{x)h' {x) 

f" (x) = -^3^^ (30) 

gi (x) p» (x) 

where, with the help of (gH) and (P7|) . 

A (x) - (x - 5: (x)) f [1 - * (s)] ds - g, (x) (1 - $ (x)) 

= g, (x) (x) (x - (x)) - (x2 - xgl (x) + g, (x)) (1 - $ (x)) (31) 
B{x)^-{x- g^ (x)) ^ $ (i) di - 5* (x) $ (x) 

= 5* (x) (x) ((?; (x) - x) - (x^ - xg'^ (x) + 5. (x)) $ (x) . (32) 
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Let d{x) — (x) [x) and n {x) = /" {x) d{x), the indicated denominator and numerator, respectively, of 
/" (x). As a; Z, both d{x) and n{x) tend to 0. If h' happens to be differentiable, then by L'Hopital's 
rule, liuix^i f" (x) ~ lim^-i.; n' (x) /d' (x). It can be shown that A' (x) = (2 — 5" (x)) [1 — $ (s)] ds and 
B' (x) = - (2 - g'^ (x)) Ji" $ (t) dt. Therefore 

n' (x) = A' (x) / $ (t) /i' [t) dt + A (x) $ (x) h' (x) + (x) / [1 - $ (s)] /i' (s) ds - B (x) [1 - $ (x)] /i' (x) 

Jl Jx 

+ [-xg^ (x) (x) + (x) 5* (x) (x)] h' (x) + (a;) (x) /i" (x) 

- (2 - .9:' (x)) / [1 - $ (s)] / $ (t) /i' (t) dt-{2- g'l (x)) / $ (t) [1 - $ (s)] /i' (s) ds 

Jx Jl Jl Jx 

+ [A (x) $ (x) - B (x) (1 - $ (x)) - (x - (x)) (x) (x)] /i' (x) + gl (x) p. (x) h" (x) 

- (2 - (--c)) gl {x) p, (x) /' (x) + • h' (x) + 5*' (x) p, (x) /i" (x) 



and so 



lim /" (x) — lim 



(2 - g': (x)) gg (x) p. (x) (x) + gg (x) p. (x) fa^' (x) 
(2g* (2:) - x) gl (x) p* (x) 



lim 



-/ (x) + lim' ^ ' 



I (x - 0^^^ on (/, 0), (x) = I |/|^/^ on [0, 00) and /i (x) = on (-00, /]. This 



x^i 2g^ (x) — X 2g^ (x) — x 

Define the function h (x) 

function is clearly Lipschitz. Note that h" (x) = on (Z, 0). We now consider the same assumptions from 

Theorem [5] and see that lim^^-^j I/' (2;)! < ^ ll^i'lloo ^^"^ linia;^/ 2g^ (x)--x ~ have thus found a Lipschitz 

function h for which limj^-j.; |/" (x)| = 00. 



Remark 10 From the above discussion we can't expect to have a universal bound on the second derivative 
of f unless the support of the target r.v. is (— cx),cx)). This is consistent with the known NP bound in 
Wiener-Poisson space developed in \30\/ . where Z was Normal and hence had (—00,00) for support. For the 
rest of this subsection, we will then assume that I = —00 and u — 00. 

Theorem 11 Assume that g* is twice differentiable and g"(x) < 2. Suppose too that is 
bounded as x —00 and as x ^ 00. Then the solution f of the Stein equation {IT^ , for a given test function 
h with ll^i'llo^ < 00, has second derivative bounded as follows: 

ll/"lloo<fc||/^'lloo (33) 

where the constant k depends on Z alone, and not on h. 

Proof. Recall the functions A and B in (I5T|) and ([5^ . Using Lemma 7 in [B] (note that $ there is defined 
as the upper probability tail), A (x) < and B (x) < 0. Therefore, from (PH)) . 



I/" (2^)1 < 



-A{x) 



gl {x) p., (x) J I 



<^>{t)dt- \\h'\ 



-B{x) 



gl (x) p., (x) 



r[i-$(s)]ds.ii/i'iu+i^ 

J X 9* 



\h'{x)\ 
(x) 



(x)p. (x)|/" (x)| 



< 



- (x - g: (x)) / [1 - $ (s)] ds + g, (x) (1 - $ (x)) 



$ {t) dt 



{x^g',{x)) / $ (t) d< + g, (x) $ (x) 



pu 

/ [l-$(s)]ds + g,2(x)p,(x) 

J X 



g4x)(l-$(x)) ^ (t) dt + g, (x) <S> (x) J [1 - $ (s)] ds + g^ (x) p^x) 

g* (x) (1 - $ (x)) (g, (x) p, (x) + x$ (x)) + g, (x) $ (x) (g, (x) p, (x) - x [1 - $ (x)]) 

+ ff* (a;) (2;) 
2ff* (a;) p* (x) . 
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Due to the continuity of and conditions of Assumption B when Z = — oo and u = od, there is some toq > 



such that g^,{x) > toq for all x e M. Then, < 



2iife;iic 

9 Ax) 



< 



\h'\ 



kWh'l 



One might think at first glance that the conditions of Theorem [TT] are too restrictive. However, a 
closer look will show that they are all satisfied by members of the Pearson family having (— (X),c») as its 
support. Examples are the Pearson Type IV, Normal, and Student's T distributions (see Table [1] to check 
the conditions). 

Assumption B' Along with Assumption B, the following hold. 
1. 5* is twice differentiable and < 2. 



2. lim^^ic 



< OO. 



4 NP bound in Wiener space 

From the results in subsection 13. II all solutions of the Stein equation belong to the set J'-u = {/ G C^(/, u) : 
1 1 /'I loo < fc}, where the constant k depends on the distance d-n used (and so it implicitly depends on the set 
H). 



Theorem 12 (NP bound) Let d-u he dw or dpM- Under Assumptions A and B, 
dn {X, Z) < kE \g, (X) -gx\< kE \{g, (X) - gxf 
< k 



E 



9* (xy 



E 



9* [zy 



E[g,{X)gx]-E[g, {Z) gz]\ + \E[gl] 



E[.9|]|. 



Let G, {x) he an antiderivative of g^ (x). Under Assumptions A' and B 
dn {X, Z)<k 



E 



9* (xy 



E 



9* (zy 



E [XG, (X)] - E [ZG, [Zy I + |E [g\] - E [.g|] | . 



(34) 
(35) 

(36) 



In hoth statements, k is a finite constant depending only on Z and on d-u ■ 

Proof. The first bound in (|34l) follows from (|21l) and Theorem [S] The second bound follows from Holder's 



Inequality. Let A 



(5* {X) ^ gxy 



1/2 



Since (5* [Z) — gz) =0 a.s. 



E 



5* [Xf] ~ 2E [g, {X) gx] + E [gl] ~ (e [5, [Zf] - 2E [.9* {Z) gz] + E 



and ([35]) follows. From and Assumption A' we have E [g*{F)gp] ^ E [FG^ (F)], which proves ■ 

The first inequality also follows from Theorem 1 and equation (19) in Kusuoka and Tudor [10^ The setup 
in their paper involves functions b and a. The function b is any function for which J^" b {x) {x) dx = 
along with a few other mild conditions: 6 > near I, b <0 near u, bp^ is continuous and bounded on {l,u). 
They then defined a{x) = 2 Jj^ b (y) (y) dy/p^ (x). Then for W a standard Brownian motion, the SDE 



dYt=biYt)dt+y/a{Yt)dWt 



(37) 



has a unique Markovian weak solution with invariant density p*. With a and b as given above, from Theorem 
1 in [10], 

a{X) 



d-H (X, Z) < kl 



{DX, DL-^ {b {X) - Eb (X)}) 



k\Eb {X)\ . 



(38) 
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If we take b {x) = —x, it follows that a (x) = 2g^, (x). If X is centered, the right-hand side of ([55)) quickly 
reduces to fcE \g^ {X) ~ gx\- 



While the results in 10 appear more general, taking b (x) — ~x suffices. A careful analysis will reveal 
that the proofs of their main results depend only on the density and the choice of b. While each choice 
of b arguably yields a different diffusion process y, the invariant density is still p*. Their analytical proofs 
are in fact independent of the stochastic differential equation (1371) and the diffusion process arising from it. 
For this paper, we only need comparisons with the law of the reference variable Z . To this end, knowing the 
density will suffice. The computations using b {x) — —x and a [x) — 2g^ (x) are much easier and this is 
reflected in the simplicity of (|34l) compared to 



Furthermore, as shown in the next theorem, the bounds we get from taking b (x) = ~x (see Theorem [T^ 
are tight. Indeed, nothing is lost by choosing b this way. 

Theorem 13 (Law Characterization) X Z if and only if all of the following are satisfied. 



1. E 



{xy 



= E 



9* {zy 



2. E [XG^ {X)] = E [ZG^ (Z)] 

3. E[gl] =Er^2- 



Proof. If the three conditions are satisfied. Theorem [T2] implies d{X,Z) = 0. 

Now suppose X Z. They then have the same density so 1 and 2 immediately follow. We next 
prove that gx ^= gz, imitating the technique Nourdin and Viens used to prove ([T2l) (see Theorem 3.1 [18]). 
Let / be a continuous function with compact support, and F any antiderivative. 



E [/ {X) gx] = E [XF (X)] = ^ [xp., {x)] F {x) dx 



-F (x) / yp^ (y) dy 



X—>1 



fix) 



VP* (y) dy 



dx 



f (x) (x) dx^E 



fx yp* iv) dy 
P* (X) 



so gx = Jx yp* {y) dy/p* {X) a.s. This has the same distribution as yp* {y) dy/p*, (Z), equal to gz a.s.. 



so 3 then follows. 



Remark 14 We see that E 



.9* [Zf = E [ZG., (Z)] = E [g|] (see Lemma\^. Thus, for X to have the 



law as Z , it is necessary and sufficient that E 



gif {X)^ , E[XG* {Xy\ and E \g\'\ (which a priori need not 



be all the same) are all equal to E \_g%\ ■ The three conditions in Theorem \13\ are stated in their current form 
due to the symmetry involved. 

That E (Z)^ = E [ZG* (Z)] = E [(]%] are all equal depends on the specific structure of Z itself, and 
it is rooted in how g^, (and thus G■^, as well) is defined in terms of the law of Z. Specifically, it is because 
(7* (Z) = gz that we are able to use the integration by parts formula (jlOp on g^, (Z). If we evaluate the 
function g^, at the random variable X , we cannot expect {X) to be equal to gx because 5* is an object 
that "belongs" to Z. However, if X and Z are to be "almost" the same in law, we would expect X to "almost" 



satisfy the same relations/equations for Z, e.g. E g* {X) 



'E [XG* {Xy\. If g* is a polynomial, then 



this amounts to checking that the moments of X satisfy the same conditions met by the moments of Z. 
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Granted, this method of moments is not sufficient. Hence, the need for condition 3, E [gx] — E , in 
Theorem [T51 



The following versions of Theorem [T^] and Theorem [T^ for sequences are useful. 



Corollary 15 X„ — > Z m distribution if all of the following are satisfied. 



1. E [.g, (X„) 

2. E[g,{Xn)gx.. 

3. E [g2.J ^ E 
Corollary 16 X„ 



E [g* (Z) g^] (under Assumption A) 
E [ZG* (Z)] (under Assumption A') 



Z in distribution if g^ {Xn) ~ gx„ — > m L^{Q,) 



Remark 17 If we normalize so that YaiX = VarZ, condition 3 in Theorem \13\ can be replaced byYavgx = 
Yax gz since E [gx] — VarX. This also allows us to replace the term |E [gx] ~ E [g^] \ Theorem by 
\Va.r gx — Yai gz\- In Corollary ] 151 we can replace condition 3 by Y&i gx„ — >■ Var g^ i/E [A"^] — > E [Z^] . 



If Z is Normal with variance so (y) — cr^ , G* (y) = a^y and 5^ = ct^. If Var AT = cr^, then 



dn {X, Z) < ky^\a^ - cr^l + a"^ |E [X^] - E [Z2]| + |Var gx - Yargzl = k^Y'Avgx 



(39) 



where fc = 4 if d-u = and fc = 1 if d-j^ = dw- This retrieves Theorem 3.3 in |15) . If we have a bound 
on Var f;x, this may be used to bound the distance. A Poincare-type inequality may be used in this regard. 
See |17| (also for an explanation of the notation used below) where they use such a bound on Var gx to get 
the following result: 



dn {X, Z) < 



fcVlO 

2cr 



E 



\D^X (E)i D^X\ 



1/2 



E 



\DX\\ 



1/2 



(40) 



This was used in [T7] and [30] to prove CLTs for functionals of Gaussian subordinated fields (applied to 
fBm and the solution of the 0-U SDE driven by fBm, for aU H e (0, 1)). 



4.1 Convergence when (7* is a polynomial 

Many of the common random variables belong to the Pearson family of distributions, all of whose members 
are characterized by their g^ being polynomials of degree at most 2, i.e. (y) = ay^ + /3y + 7 in the support 
of Z. Some member distributions in this family are Normal (g* is constant), Gamma (17* has degree 1), 
Beta (17* is quadratic with positive discriminant). Student's T-distribution (g* is quadratic with negative 
discriminant) and Inverse Gamma (g* is quadratic with zero discriminant). 

Refer to [B] and [55] for more information about Pearson distributions, and [5] for Stein's method applied 
to comparisons of probability tails with a Pearson Z. From Remark 01 if the support of Z is unbounded and 
g^ is a polynomial, then Z is necessarily Pearson. If Z has bounded support and g^ is a polynomial, may 
have degree exceeding 2 and in this case, Z is not Pearson. 

Corollary 18 // 5* is a polynomial .g* (x) — Y^^k^Q CLkx'' , for the convergence Xn Z in distribution, 
conditions 1 and 2 in Corollarv \15\ can be replaced by these conditions (respectively): E [X^] — > E [Z*^] for 
k — 1, . . . , 2m, and E [Ar,^gx„] — > E [Z^gz] for k = 1, . . . , to. Under assumption A' , the two conditions can 
be replaced by E [x!^] E [Z'^] /or fc = 1, . . . , max {2m, m + 2}. 
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Proof, g'^ (x) has order 2m while xG■^, (x) has order m + 2. The matching moments ensure condition 1 in 
CoroUary (TS] is satisfied, and under Assumption A' also condition 2 is fulfilled. ■ 



Suppose (x) = X^fcLo ^kx'^- Note that 



E 



9* {zy 



E 



\k=0 



fe=0 \j=0 



while 



E 



k=0 



We noted earlier that E 



(Z) and E [ZG* (Z)] are equal. While the polynomial coefficients of the 
different moments of Z are different, and more moments may be involved in one expression compared to the 
other, the coefficients and the moments themselves should take care of this apparent difference to ensure 
equality under the expectation. 

Suppose Z is Pearson with gz = g* (Z) = aZ^ + I3Z + 7. Using Lemma [Sj we can prove the following 
recursive formula for the moments of Z: E [Z'+^j = j^E [Z''] + j^E [Z"-^] . Therefore, 



E [gz] - E [Z^] 
2E [Zgz] = E [Z^] 
3E [Z^gz] = E [Z^] 



1 - a 



2/37 



(f - a) (1 - 2a) 

6/327 + (1- 2a) 372 
(l-a)(l-2a)(l-3a) 



and 



IE [.g|] 



P^-f (1 ~ a) + j^l - 2af 
(f - a) (1 - 2a) (1 - 3a) 



Var,. = E [gl (Z)] - (E [9. (Z)])^ = /^^T d - a) V 2a V (1 - 2a) ^ 
^ ^ (l-2a)(l-3a)(l-a)' 



(41) 
(42) 



Corollary 19 Suppose Z is a Pearson random variable and for the sequence {Xn}, VarX„ = E[X,'^J 
IE[<7Jf„] ^ following are sufficient conditions so that Xn Z in distribution. 

1. When Z is Normal (a = [3 = Q), Ysngx,, — !• 0. 

2. When Z is Gamma (a = 0), Yargx„ — 13'^ J and 

• under Assumption A, E [X„gx„] ■ 

• under Assumption A' , 2E [XngxJ = E [-'^,1] ^ 2/37. 

3. In the general case where a 7^ 0, Varg^ — > ^ 7/"'"o ".li"*" o" Ji <ind 

^ ' ' (1 — 2q)(1 — 3q)(1 — a) 

• under Assumption A, 



2E [XngxJ,E[X^] 



2/37 



• under Assumption A' , 
2E [X„gxJ = E K] - 



(l-a)(l-2Q) ■' 



2/37 
(1-q)(1-2q;) ' 



id 3E[X25^J ,E[X4] ^ 
and 3E [X^g^j^J = E [X^] 



6;3'^7 + (l-2a)37^ 
(l-Q)(l-2a)(l-3a) ' 



6;3^7+(l-2a)37^ 
(l-a)(l-2Q)(l-3Q;) ■ 
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Proof. Apply Corollary [T51 directly. ■ 

The first statement is the yersion for sequences of Corollary 3.4 in [18^. Alternatiyely, we could replace 
Ysxgx^ — > by E [g^,,] ^ 7^- -for the Gamma convergence, we can replace Ya,rgx„ — >■ /3^7 by E [.gx„] ^ 
/3^7 + 7^. When a 7^ 0, we can work with (1411) instead of p2)) so the statement will be in terms of 

Iff r„2 1 , /3^7(l-a)+7"(l-2a)^ 
LyX„J ^ (1-q)(1-2q)(1-3q) ■ 

The next result follows from Corollary [TBI 

Corollary 20 Suppose Z is a Pearson random variable. Xn Z in distribution if gx„ ~ ctX^i " PXn ~^ 7 
in (fl). 

4.2 Convergence in a fixed Wiener chaos 

When X is inside a fixed Wiener chaos so X — Ig (/), we have more structure ayailable. For example, 
{DX,-DL-^X)^ = i \\DX\\l. Therefore, if Z ^= J\f{0,<T'^) and E [(/, = a^, ^ gives us the 

bound 



d-H {X, Z) < k^Y&igx < kJvsiT Q \\DX\ 



One may then use bounds like 



Var ( i \\DX\\l 



(a) 1 



\DX\\l-qa' 



(43) 



to further cap the distance. Equality (a) follows from E 



E [gx] = and inequality (b) from 



Lemma 3.5 in [15]. These are quite important and known results which yield CLTs for functionals on a fixed 



Wiener chaos. For instance, if we have a sequence {Xn} — {Iq (/«)} where E {fn)Y 
following conditions are equivalent : 

1. Xn — > Z in distribution; 

2. E [X*] ^ 3a4; 

3. Il/n <8)r /n|lH»(29-2r) ^ for aU T = 1, . . . , (J - 1; 



(7^, then the 



4. \\DXn\r^ ^ qa^ in L^n); 

5. ||i?2x„®ii?2^„||J^^, ^OinL2(r!). 



See [H] for (1) <^ (2) ^ (3), [lOj for (1) <^ (4), and [H] for (1) ^ (5). These in some sense 
highlight the tightness of inequality (p6l) with the help of bounds like pO|) and ([43]) . 



Corollary 21 IJ Xn — Iq (fn) with q > I, then condition 3 in Corollarv \15\ can be replaced 6t/E ||Z3X„| 

q'^ [gUz)]. 

Proof. This is a direct consequence of {DXn, -DL-'^Xn)^ = \ \\DXn\\\ and E = E [gl {Z)] . m 

From this and Corollary [201 we have the following result for the convergence in a fixed Wiener chaos to 
a Pearson random variable. 
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Corollary 22 Let Z be Pearson with (z) — az^ + /Sz + j in its support. Fix q > 2. Suppose Xn = Iq (fn) 
and E [X^] — > If — qaX^ — qfiXn qj in L^ (51), then X„ Z in distribution. 

Remark 23 Special cases of the above corollary are known results. 

• Let Z be Normal with variance 1, i.e. g^ {z) — 1. Suppose E [X"^ — >■ 1- Then Xn Z in distribution 
if \\DXn\\% -^qinL^ (17). See l2Ql 

• Let Z be Gamma with g*{z) = (2z + 2v)^, i.e. /3 — 2 and 7 — 2v, where the parameters are chosen 
for consistency with the discussion in '161. Suppose E [^^] — > 2v. Then Xn Z in distribution if 
WDXnWl - 2qXn ^ 2qv m L^ (n). 

The result in the first item of this remark is known as the Nualart— Ortiz-Latorre criterion. In [551, the 
authors used it to prove that 

C%/7Vln(iV) (Hn - h) > 7V(0, 1) 

where is an estimator of the Hurst parameter H for fBm when He (51 -j) (see [55] for details). 

4.3 Bilinear functionals of Gaussian subordinated fields 

Let Xt be a centered Gaussian stationary process with covariance function C (t) =K [XoXt] = E [XsXg+t]. 
Let / : M ^ R be a non-constant function such that with Z ^ 7V(0, C(0)), E[|/(Z)|] < 00 and 
E I/" (Z)l** < 00. Write ^/ = E [/ (Z)]. It was proved in [17 and [30^ (under mild conditions) that 

'-^ ^ / (^^) - — ' -^(0, S2) 

where V{T) is a normalization function with specific properties. In particular, for the case of the increments 
of fBm {Xt = B/^i - Bf) with Hurst parameter H £ (1/2, 1), we have V{T) = T" and = 2(E[Z/(Z)])^ 
As an ilustration of how to employ this tool, and to emphasize the advantage of using inequality (j35p over 
((M)) . we will prove that Ft :— (-^)^ — E (■^)^ converges to a (centered) chi-squared r.v. as T — >■ 00, for 
the case of the increments of fBm when H e (1/2, 1). 

Preliminary computations and notation: 

Following the setup in 19 (Section 5.1.3), we can write B/^i - = /i(i^ij(t, •)) where Knit.s) = 

CRsi^^ J^^^{u — s)^^2u^^2dul[Q t+i] and ch = i3(2-2H~h-^) " "^^^ integral /i is with respect 
to a Wiener process W generating the same filtration as , with the two processes related by Bt — 
jlKH{t,s)dWs. 

Let's make the simplifying assumption C (0) = 1 so Z, Xg ^ A/" (0,1). To simplify notation, we will 
write /j,, for /j„ and Kt for KH{t,-)- Also let Sj = L^{[0,T]). Then DXt = Kt and E[A^At] = 
E[/i(Ks)/i(Kt)] = {Ks,Kt)^^ ^Ci\s-t\) Cst. Define for T > and s,t £ [0, T] the functionals 

F,t = {f{x.)-^lf){f{Xt)~^if) 

Ft = (^^) = J^^ if (Xs) - lif) if (Xt) - ^Lf) dsdt = j^^ Fstdsdt 

Ft — Ft — E Ft ■ 
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Let s = (si, . . . , sp) G [0, T]^ . We will use e(si, Sj) to denote a nonnegative integer exponent indexed by a 
pair of variables from s. Define 



ds = 



1 



n c^s- 



where inside the integral is the product of Q = (^) factors. For example, if P = 4 and s = (s, t, u, w), then 

1 



LiT) 



^e{s,t) ^t(s,u) /-it(s,v) /^<i{t,u) ^t(t,v) ^c(u,v) 
^st ^ su ^sv ^tu ^tv ^uv 



dsdtdudv. 



(44) 



Recall that C^^ ' = (C(|s — t|)) ' , so the integration in (H^ is being done only on the subscripts and 
not on the superscripts (since these are fixed exponents indexed only by the variables over which we're 
integrating) . 

Proposition 24 With the previous notation, 

1. For q>l, take Cqq\ = E [Hq (Z) f (Z)] where Hq is the q^^ Hermite polynomial. Then, 

oo 

f{X,)^f,,+Y^CqIq{Kf^). 

9=1 

2. Fst has Wiener chaos decomposition 

oo n— 1 oo n oo 



CrnCn-mIn (xf ™ ^ xf ("-'"^) + ^ ^ ^ d (m + r, n - m + r, r) ( 



K?™(8)K 



n—2 m—1 

oo n oo 



n— ni—0 r—1 



= e (to, n, r) d (to + r, n — m + r, r) !„ ^ 

n— m— r— 

w/iere d (fc, j, r) = CfeC^r! (^) (^) , and 

e (m, n, r) - 



K?"(8)K 



if r = and to G {0, n} 
f otherwise 



3. For Q < a, b < n. 



5(n— a) 



(„-b)\ 



where the summation is taken over all p for which max (0, a + 6 — n) < p < min (a, &) 
4- Fix an integer P > 2. Let S = e {si, Sj) be the sum of the exponents in L (T). 

• If S > P/2, then limT^oo L (T) = 0. 

• If S — P/2, then limT->oo L (T) < oo. 

Proof. To prove the first point we expand / in terms of Hermite polynomials: 

oo 

f{z)=C^+Y,CqHq{z) 
9=1 



f-ip /^a—pf^b—p^n—a—b+p 
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with Cqq\ = E[i7g [Z] f [Z]] for all q > Q. Since for any h e L'^{\Q,T]) we have the relation Hq{h{h)) = 
Iq{h'^'^), then the result follows. 

For the second point we have that Kf'=(^,Kf = (K*, K,);^ (^Kf Kf^'"'^) = (xf ^'^"''^ Kf 
Therefore, from the previous point and the product formula (j2|. 

oo oo oo 

fe=i j=i 



fej = l r=0 \ / \ / kj=l r=0 



Write c. 



(fc,j,r) = d(fc,j,r)CJ,/fe+,_2. (k?('=-'-) 0Kf(^"-'-)) where d (A:, j, r) = c,c,r!(^)pj. Then 

OO OO kAj oo oo oo oo kAj 

= E E E ^) = E E 0) + E E E ^) • 

fe=lj=lr=0 k=l j=l k=l j=l r=l 

Applying Fubini's theorem for sums we have, 

C30 oo OO OO OO OO 



k—1 j — 1 m—1 j—1 m—1 n—m-\~l 

oo 71— 1 



n—2 m—1 



and 



Therefore, 



oo kAj k j oo k k k k oo k oo 

EE-EE+ E E-EE+E E =EE- 

j—1 r—l J — 1 r— 1 j—k--\-l r—1 r—1 j—r r—1 j—k-'rl r—1 j—r 



oo oo kAj oo k oo oo oo oo 

^^^cst = XIEE"^"* (^'■^''^^ = E E E ^ 

fc— 1 J — 1 r—1 A:— 1 r—1 j—r r—1 m—0 n—m 

oo oo n 

= ^ ^ ^ Csi (m + r, n — m + r, r) 

r—1 n— m—0 

oo n oo 

(n—m) 



^ ^ ^ d (m + r, n - ,71 + r, r) CIJ^ (xf ™ Kf 

n— m—0 r—1 



establishing the second point. 



Point 3 requires counting the possible combinations in the inner product. Note first that the symmetric 
tensor product Kf'iKf ^"""^ has (^) distinct terms. Take any particular term a and list down all its n 
factors Kt and Kg in the order in which they appear. Now take any term j3 from K^''(g)K?^" and list 
down all its factors (in order) below those of a. Let p be the number of (Kt, K„) pairings. Thus, the number 
of pairings of the type (Kj, K.y), (K.;, K„) and (Ks, Ky) are a — p, b — p and n — a — b + p, respectively. 
Finally, the number of pairs (a, /?) which have p matching Kj and K„ is b-pn-a-b+p) ■ 

Finally, to prove the fourth point we make use of Proposition 3 (point 4) in 'W' which states that, 

/ rp2H 

ClAdsdt^O 

T2 



" I y(2-2ff)(£-l) 
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so 



lie: 









i/Q 








i/Q / 


\I 




ds 




j.P-2 f 


dsdt 


= o( 















st\\Q 

By the generalized Holder inequality for integrals, 

1 



' rp(P-2+2H)/Q 
Xi2-2H)iQe-l)/Q 



L{T) 



o 



and the result follows. 



n Cs.s 



^PH n 



(Si,Sj) 



< 



PH n ll*^^ 



rpPH 



j^(P-2+2H)/Q 



TPH 11 rp{2-2H)(Qe(s,,Sj)-l)/Q 



■e(si,Sj) 



o 



1 



y(l-H)(2S-P) 



Theorem 25 Suppose / 



is o/ dass such that 2(E[Z/(Z)]) ^ 0. T/ien as T ^- oo, 

-r-, Law 9 
i^T ^ X ■ 



Proof. Throughout this proof we will use the symbol ~ to relate two expressions which have the same limit 
as T — ^ oo. Also, k will represent a constant whose value may change from one equation to another and 
might depend on the summation indices. 

Since is a- member of the Pearson family, we can make use of Corollary [1^] to prove this theorem. 
This will require proving convergence of the second moment of Ft and of g p^. , plus the moment convergence 
of FtQEt- We take note of the following facts: Y? = 2c? (ci E[Z/(Z)]), and T-"^" jj,^ Ctsdtds ^ 2 as 
T — > cx) (see Proposition 3 in [30]). 

• Convergence of the Second Moment of Ft'- 

From point 2 of Proposition [Ml 



oo n oo 



DFst = ^ ^ ^ ne (m, n, r) d (m + r, n - m + r, r) /„_i (Kf"' (g) K' 



(n-m)\ 



and 



71—1 m— r— 



AT 



(45) 



DL-^Fu, = I] II I] e (Af, N, R) d{M + R,N~ M + R, R) In-i (k^*^ (g> K?(^-^^)) Cf, 



Ar=l M=o R=0 



(46) 



Each indicated multiple integral of order n — 1, with kernel in i^**", is to be interpreted as 



lu-i k: 



0m 



K 



® (n— m) 



TO 

n 



/„_i (Kf ® 

l^{m,n,a) 



E 



K 



5(n— 1 — a)\ n(m,a) 



(47) 



where fc(™--)fir) = ^ 7^^^ if a = m - 1 > 

' [n — to) Kt it a = m 
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Therefore, E 



(DFsu-DL-^F^,)^ 



will be a summation containing terms of the form 

{N-l^A)\ r-{M^A)\ 
j uv / 

o(m,a) n{M,A)\ /-ir /-iR 
st 1 uv I ^ sV^ui 
I S 



t I '■St 



I J^»(ti-1) 



where from to ([50)) . we used the third point of Proposition [ 



^ uv 



(48) 
(49) 
(50) 



Observe that ( f ^ , 



^(m,a) r{M,A)\ 



is any of Cs„, C^^, Ct„ or Ct„. In ([50]) . the exponents of all six types 



•■St I '■I't' / 

of correlations then add up to S = n + r + R — N + r + R (at this point, always take N = n, 
otherwise the term is 0). In the summations P5|) and appearing in (^DFst, —DL^^Fuv) ^, since 
e (m, n, r) e (A/, A^, i?) = if 5 < 1, the remaining terms are those for which 5 > 2. Therefore, using 
the moments formula of Proposition [SJ 



E [F|] = E [gp^] = E {DFt, -DL-^Ft) 



E 



{DFsu-DL-^Fuy)^ dsdtdudv (51) 



S4^^ ^4 ^ ^^^"^ ® Kt) , d (1, 1, 0) h (K„ ® K,))^,] dsdtdudv (52) 



where in (|52p . we applied point 4 of Proposition [Ml (P = 4) on (|50p and (|?T|) : for those terms 



contributed by E 



{DFsu-DL-^Fu,)^ 



where S > 2, the limit is 0. The limit in (I5ip is the nonzero 



value we get for the remaining case 5 = 2; specifically, (n, iV, m, M, r, i?) = (2, 2, 1, 1, 0, 0). Since 

h (K,) /i (KJ Ct, + /i (Kt) h (K„) 



(/i(K,®KO,/i(K„®K,))_, 



Ji (K«) /i (K„) Ct, + h (Kt) /i (K„) a„ 



(53) 



then 



Therefore, 



[^t] - j^^ \ [^CsuCtv + 2CsvCtu] dsdtdudv ^ ^ (2)' = 2. 



E i^^ ^2 



Convergence of FtQEt' 

Notice that 

E [Frgp^] = E 



Ft9Fi 



E 



FtQFt 



since 



E 



Ft 



/ E [Fst] dsdt = y '^^I'YJ [ c:,dsdt ^ / c. 



St 



dsdt 



2c| 



Now we need to investigate 
E 







Ft 9 Ft 


.We 



{DF,t,-DL-^F,,,) F^, 



dsdtdudvdwdx. 



(54) 
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The expression inside the expectation is a summation with generic term 



>(n— m) 



(55) 



where a G {m - 1, m}, A e {Af - 1, M} and (4 , fi 
(we used (H71) here). 



(m,a) P^A^.'^)^^ ^^^j^ ^^^^ ^^^^ 

Si 



(nA7V)-l 

E 



K 



5K. 



5 ^n' — m'^ 



l{n'^n+7V-2-2/} 



5: ^ 

{2/=n+W-n'-2} {eSA/} 



E\ ^ (-^c(s,u) /-it(s,v) (-i<^{i,'ii) /^i^{t,v) ^(^(s.w) /-i^(s,x) /-<i^^ 

/ J ^su ^ sv ^tu ^tv ^ sw ^ sx ^tw ^tx ^uw ^ux ^vw ^vx 

{2f=n+N-n'~2} {eeAf,ieBf} 

(56) 

where Af is the coUection of exponents e{t,u) = p, e(t, w) = Si — p, e(s,u) — S2 ~ p and e{s,v) = 
f — Si — S2+P, withp, 5*1, 52 running over all integers such that max(0, + 6*2 — /) <p< min(S'i, iS'2), 
max(0, n— l — a— f)<Si<a and max(0; N — 1 — A — f) < Si < A. is defined similarly, just 
applying recursively point 3 of Proposition 1241 so for instance we have £,{t,w) + £,{t,x) + ^(s, w) + 
^(s, x) + £,{u, w) + ^{u, x) + ^{v, w) + ^{v, x) = n' . From this point on, we take 2f = n + N — n' — 2. 
Combining ([55t and (|55|) . we see that the expectation inside the integral of (|5^ is a summation 

oonoooo A''c300o n' oo 

EEEEEEEEEE E E 

n=l m=0 r=0 N=l M=0 R=Q n'=0 m'=0 r'=0 a, A {2f=n+N-n' -2} {eeAf.^e Bf} 

with each term consisting of 15 different types of correlations whose exponents add up to S* = / + n' + 
1 + r + R + r' . From this it follows that n + N + n' + 2 {r + R + r') = 2S. Systematically listing down all 
possible values of the indices in this summation will show that e (m, n, r) e {M, N, R) e (m', n' , r') ~ 
when 5 < 2. For those terms in the summation for which S* > 3, the limit they contribute in (I54p 
is (using the fourth point in Proposition 1241 with P ~ 6 variables s, i, u, v, w, x). For those terms 
having 5* = 3, a careful consideration of the indices such that e (m, n, r) e (M, N, R) e (m', n' , r') > 
leads to either (m, n, r, M, N, R, m' , n' , r') = (1, 2, 0, 1, 2, 0, 0, 0, 1) or (1, 2, 0, 1, 2, 0, 1, 2, 0). Therefore, 
we can continue (l54l) as 



E 



FrgPi 



2cl 



= L1+L2. 



E [(/i (Ks ® Kt) , h (K„ ® K„))^^] C'wxdsdtdudvdwdx 
/ E [(/i (K^ (g) Kt) , /i (K„ (g) K^,)) , /2 (K„ (g) K^c)] dsdtdudvdwdx 
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Using ([55)1 . we have 
2cf 1 



T-^H J^^ 



J ^ Cyj^dwdx ■ J ^ E [(/i (K, (g) Kt) , /i (K„ (g) K„)) dsdtdudv 



2c? 



rp4 



^ dsdtdwdi; ^ ^ • (2) • (2)' = 2 ( 91 



1 /" 2CsuCtv + 2Csj,Ct, 



2c 



S2 



,2\ 3 



To compute L2, we use ((53)) again so 
E [(/i (K, <E> Kt) , h (K„ ® K^))^ I2 (K^ ® K,)] 

= ifi [/i (K,) h (K„) /2 (K^ ® K,)] Ct, + [h (Kt) /i (K^ (K^ ® K,)] 

+ ^E [h (K,) /i (K„) /2 (K„ ® K,)] Ct„ + ^E [h (Kt) /i (K„) h (K^ ® K,)] 
The first expectation simplifies to 



r=0 



E [h (K,) /i (KJ /2 (K„ ® K,)] = ^ r! ( ^ ) ( M E [/2-2r (K, ®, K„) /2 (K,„ ® K,) 



2! ^Ks(g)K„, KujgjKjj^^j — 2 — CgwCux + C's^C'-u 



so 



2cf 



2cf 

^6 



/ E [(/i (K^ Kt) , h (K-u (g) Ki,)) - /2 (K„ (g K^)] dsdtdudvdwdx 



{CswCux CsxCuw) Ctv ~t~ {CtwCux ~t~ CtxCuw^ Csv 
{C swCvx CsxCvw^ Ctu ~t~ {CtwCvx ~t~ CtxCvw^ C s' 

„6 



dsdtdudvdwdx 



2I]6T6^ 



/" 4c 3 

• 8 y CswCuxCtvdsdtdudvdwdx -> ^ (2) = 4. 



Finahy, 



; [FtQFt] ^ (2 + 4) - 2 = 4. 



Convergence of the Second Moment of Qft- 



E 



Ji8 



X (Kr' «Kf (k^^' «Kf 

= ^ (Kf^ ® Kf ("-^-''') /at^i (k^ ® K?(^-i-^)) [k^' ^ k! 



(57) 



{DFst,-DL-'^Fuv)^^{DF^^,-DL-'^Fy^)^^ dsdtdudvdwdxdydz. (58) 
Inside the expectation, a generic term is 

(Kf " ® Kf ("-™)) Clt, /^-i (k?*' ® K?(^-^^)) C« )^ 



X /iv'-i k: 



6DKf ("'-^-^')) (fit"-\f(f-^))^^ (fir' '^'Ifif ' ^'))^ 



ryR r^r /^iR 
^ st^uv^wx^yz ■ 
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The expectation of the product of the four multiple integrals is of the form 
E [/„_i (a) /jv_i (6) /„'_! (c) In'~i {d)] = ^>^ {a(§pb, c(§qd)^^( 



E 



(^DFgt,~DL ^Fuv) {DF^^j., ~DL ^Fy^)^^ then is a summation 



oonooooA^oooon'oooo A^' oo 

EEEEEEEEEE E E E E 

n=l m=0 r=0 Ar=l A/=0 fl=0 n' = l m'=0 r'=0 Ar' = l Af'=0 i?'=0 a,A,a',A' {n+N-2p=n' +N' -2q} 

consisting of (2) — 28 different types of correlations whose exponents add up to S* = p + q + 
{n + N -2-2p) + 2 + {r + R + r' + R'). Along with the condition n + N - 2p = n' + N' - 2q, 
we have 25 = n + N + n' + N' + 2{r + R + r' + R'). By a careful consideration of the indices, 
e{m,n,r)e{M,N,R)e{m',n',r')e{M',N',R') = if 5 < 3. By point 4 in Proposition [Ml using 
P = 8 for the number of variables in ((58)) . we see that the integrand terms for which S > A contribute 
nothing to the limit as T ^ 00. If 5 = 4 and e (m, n, r) e (M, N, R) e (m', n' , r') e {M',N',R') = 1, we 
only have (m, n, r) = (M, TV, R) = {m' , n' , r') = (M', N' , R') = (1, 2, 0). Therefore, 

^ [sj^t] ^ ^8y8g j ^ [(-^1 K*) ' ^1 ^v))^-, {h (K^o ® K:,) , h (Ky ® K^))^] dsdtdudvdwdxdydz 

We use ([55)1 on the integrand: 

E [(Ji (K, (g) Kt) , /i (K„ <g> K,)),^ (Ji (K^ K,) , h {Ky (g) K,))^] 

.1e 
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{Ji (K,) h (K„) a„ + /i (Kt) /i (K„) + h (K,) Ii (K,) a„ + h (Kt) h (K„) C^u} 
X {/i (K^) /i (Ky) + /i (K,) h (Ky) C^, + /i (K^) h (K,) C^,^ + h (K,) /i (K,) J 



Therefore, 

^ ^^^t] ^ ^^Lf 'Ye'^^Jrs^ ^^'^ ^"^"^ ^"^"^ ^^^^^ CtyC^.dsdtdudvdwdxdydz 

~ 1^ [^1 (Ks) /i (K„) /i (K^) h (Ky)] dsdudwdy ■ ^ J^^ Ct,C^,dtdvdxdz 

We have 

1 1 

E [h (K,) /i (K,) h (K„) /i (Kj,)] =Y.J2^ [^2-2p (K, ®p K„) h-2q (K,„ (gp Kj,)] 

= E [/2 (K, (S> K„) /2 (K„ Ky)] + E [lo (K, ®i K„) /„ (K^ (»i Ky)] 

_^/ (K, 0K„) + (K„0K,) (K^(8)K^) + (Ky0K^) \ 

\ ^ ^ I j5®2 

^s?iJ^uy ^uuJ^sy ^sy^uw ^uy^sw ^ ^ 

— Csxo^Xiy ^ C^^Csy ^ ^su^wy 

and so 

A„S IP If 4c^ 

^ [Sf^] ^ t^'J^'^J^, Cs^oCuydsdudwdy ■ ^ j^^ CtyC^.dtdvdxdz ^ ^ • 3 (2)' (2)' = 12 

proving that 

E [gl^] ^ 12. 
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Let 13 — 2 and 7 = 2. We've shown that 



2-7 
4-/37 

12 = /327 + 72. 



Therefore, -Ft converges to a (centered) Gamma random variable. Specifically, since /3 = 2, it converges to 
a (centered) Chi-squared random variable with one degree of freedom (see Table [Ij . ■ 

Notice that this example is not trivial since is not a process in a fixed Wiener chaos. Also, it shows the 
importance of equation ([M)) in Theorem [T^l since equation ((M)) is very intractable in this case. 



5 NP bound in Wiener-Poisson space 

In Wiener-Poisson space, if we repeat the process before equation ((2T|) and use ([8]), the correct integration 
by parts formula, we get 



dH{X,Z)< sup 



E[/' (X)(5, (X)-.gx)]+E 



f"{F + xu)x{D^F - u)du, -DL-^X 



(59) 



It becomes evident that we need to find universal bounds on the first and second derivatives of /. Recall 
from subsection 13.21 that we only have such bounds when I — ~oo and u — oo. With this in mind, we have 
'^n = {/ G : /' is Lipschitz, ||/'||oo < fci, ||/"||oo < ^2}, where ki and k2 depend only on the distance d-u- 
The following is a generalization of Theorem 2 in [SOj (where Z was standard Normal) and an extension of 
Theorem [T2l to Wiener-Poisson space. 

Theorem 26 (NP hound) Let d-u be dw or dpM- Under Assumptions A and B' , 
dn {X,Z) <k(E \g, {X) - gx\ +E[(jx {DXf\ , \-DL-^X\ 
where k is a finite constant depending only on Z and on d-^ . 

Proof. This follows immediately from ([5^ since |(a, &)| < (|a|, \b\) and |//| < / |/|. ■ 

This upper bound was first developed for Poisson space in [22] . where was used to prove several CLTs 
for Poisson functionals. In [30 it was used to prove CLTs for Wiener-Poisson functionals. 

Corollary 27 X„ Z in distribution if both statements are true. 

1. g,{Xn)-gx„^0 mL^n). 



2. 



X {DX„) 



, \-DL-^Xr, 



in L\n). 



For convergence results inside a fixed Wiener chaos, the following preliminary computations are needed. 

Proposition 28 Let Xn — Iq (fn), with E [X^] = ql ||/n|||,®q ^ 1- Assume that for r — 0, . . . ,q — 1 and 

s = 0, . . . ,q - r, Wfn ®r /„||^»(2,-2.-=) l{s=o,r^o}u{s^o,r=o} ^ 0. Then as oo, 



1. E 



\DXn\ 
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2. ||DX„||2 zni2(f2); 



3. E 



I E [Xt] ^ 3. 

Proof. Since D^Xn — qlq-i {fn {z, ■)), we can apply the product formula (jH) to get 



\\DXn\\;^ = {DX,,,DX,,)^ = q' I (/„ (z, •)) 1,-1 {fn {z, •)) dfi (z) 



q—1 q—l—r 

E E ' 

r=0 s=0 
9 9-p 



\ 2 / \ 

q — 1\ / q — 1 — r 



l2q-2-2r-s {fn {z, ') 0^ /" ■))df^ (z) 



9'EE(P-l)'^'('^_l) ^) I h,-2,-s{fn{zr)®Ufn{zr))dn{z) 



q q—p / \ 2 

q q-r 

— rr!s! 

r=l s=0 

Also by orthogonality of chaoses, 



> 2 / \ 2 

9\ 



hq-2p-s I fn (2, •) fXlp-i /„ (z, •) d// (z) 



l2q-2r~s {fn K fn) ■ 



E 



g q— r q—B. 

= E E E 

9 q-rq-R. 

- E E E 



r/ Vi? 



q — r 



q- R 
S 



E [/2q-2r-s [fn fn) ^2g-2fl-S 



•.H=i s=o S=0 



g — r\ fq — R 



X l{2r+s=2fl+s} (2g- 2r- s)! ||/„ ®l /„||^»(2,-2.-,) ||/„ /« 1 1^,»(2,- 



2r-s) 



q-1 q~r / \ 4 / \ 4 



(r) C/) (2'?-2r~s)!||/„®^/„||^«(2,-2,-., 



In (*), we used \\g\\s^ < llslU for nonsymmetric g (this follows by a simple application of the triangle 
inequality), and Holder's inequality in the following: 



E [/2q-2r-s {fn ®r fn) I2q-2R~S {fn ®% fn)] = l{2r+s=2fl+S} (2? - 2r - s)\ {^/n®^/" , /n«)fi/n ) 



< l{2r+s=2fl+s} (2g - 2r - s)! 



fn®rfn 



j^»(2<j-2,--s) 



^»(2,-2r-a) 

fn®Iifn 



^®(2<,-2r-«) 



E 



then follows from the assumptions on the kernels' contractions, proving the first point. 



On the other hand, 



E 



\DXr,\\'-q 



= E 



\DXX,-M\DXn\\i+q^ 



= E 



\DXn\ 



-2q-qK + g2 ^ 



so ||DX„|L — > q in (fi) proving the second point. 
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For the third point we have, 



r=0 s=0 ^ / \ ' / 

q-1 q-1 q-l-rq-l-R , i\2/ i\2/ , ^2/ , nN2 



S 



E 



(i?,X„)^ (z) 



X l2q-2-2r-s{fn {z, •) /" (^^ ' )) ^2g-2-2/?-S (/« {z , •) /„ (z, •)) 

q-1 q-1 q-l-r q-l-B. , -,\2/ i\2/ , \2/ , nx2 

r=0 7?,=0 s=0 S=0 ^ ' ^ ' ^ ' ^ 

X / E \l2q-2-2r-s {xf„ {z, ■) (g)^ /„ (z, •)) l2q-2-2R-S (xfn (z, •) ®| /„ (z, •))] (2) • 



The expectation, when 2r + s = 2_R + S", is bounded by 
{2q~2r~ s- 2)! (a;/„ (z, •) «)^/„ (z, •) , x/„ (z, •) ®%fr, (z, •) 



^»(2<j-2H-S-2) 



< (2g-2r-s-2)! a;/„ (z, (z, •) 



(z, O'X'i?,/™ (z, •) 



j5|®(2g-2r-s-2) 

Modulo the constant factor (2g — 2r — s — 2)!, the integral of the expectation is bounded by 



^®(2g-2K-S-2) 



xfn (Z, ■)®rfn (z, •) 



^®(2,-2r-s-2) 



xU{z,-)®nfn{z,-) 

xfn (z, •) ^i^/n •) 



(z) 



^»(2,-2r-s-2) 



a;/n(z,-)«)i?,/n(z,-) 



j5|®(2g-2R-S-2) / ^ 



< {z, •) /n|!ij«(2,-2.-.-2) X ||||x/„ (Z, •) /„ (z, •) | |jo»(2,-2R^S~2) I 



We'll work out the first factor: 



1 1 2 / 2 

k/n (z, •) ®l fn (z, •)IU»(2,-2.- = -2) |L = / \\xfn (z, •) (g)^ /„ (z, •) |1 ^»(2,-2.-s-2) C?/i (z) 



I (/n ^r^"'" /n) (z, ') ||ij»(2,-2r-a-2) "^M (^) — 



s+1 ^ ||2 

g-2i--a-2) ""A^ || J " J " 1 1 ^ 8 (2ij - 2i-- a - 1 ) 



Finally, 



E 



{D.Xnf dfi (z) 



q—1 q—l — 7'q—l — R 
r,i?=0 s=0 S=0 



X l{2r+s=2K+5} (2g - 2r - S - 2)! ||/„ /n||^»(2,-2.-,-l) ||/« ^r""' /, 











'n ( 












r 





2 / ^ ^ 2 

s 



s + 1 r ||2 

fi Jn||^»(2g-2R-S- 



i? 



< - 1 



T- 1 

|2 



= 9' E EE^!^!(^"i)!(^~i)! 

r,fl=0 t=l T=l 

X l{2r+t=2_R+T} (2? — 2r — t — 1)! ||/„ (g)* /n||^8(2,-2,— t) X ||/„ (8);^ fn\\^tg(2q-2R-T) 

The third point then follows. 

Finally, for the fourth point, we use the moments formula Proposition [5] to get 



[X^]^-E[x'J\DXn\\l 



X {DXnf , Xn + B.xDXn) 1 = -C/„ + - (K + VK„) 
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where 



W„ =E 



x{DX,,)\X^ 



x{DX,,)\0,xDX^ 



E 



x{DX,,)\X,,{DX^) 



E 



^x'^{D^Xn) d^i{z) 



and Un — K 



xl WDX^wl 



It is sufficient to prove tliat C/„ — )• g, y„ — )■ and — >■ as n — )■ oo. 



To compute Un, note that X^ = J2l=o Z^Lo -^2g-2r-s (/« <8)r /«)• Using our expression for 



\\DX,i\\^ above, 

q q q—r q— R 

r^O _R^1 s^O 5^0 
g q q—rq—B. 

= E E E E 

r^O R^l s^O 5^0 
t/ q— r q—R 

= j2 E E 

r=o,n=i s=0 S=0 



2/\2/ \2/ 7-,\ 2 

9^ (q-R 
R 

2/\2/ \2/ nx2 

f q — r\ ( q — R. 



E [/29-2r-s (/« ®' /«) hq-2R^S (/n «)|^ /„)] 



r/ Vi? 



5 



q — r\~ f q — R 
S 



i{2r+s=2R+S} - 2r - s)! ifn®rfn, /n«)_R./, 



l{2r+s=2i?^+S} (2g - 2r - s)! {fn®rfnJn®Rfn 



i5»(2g- 



r=l s=0 



{2q~2r~s)l 



^»(2<j-2,--s) 



II/. 



2 



We can again apply Holder's inequality on the inner product, and conclude that all the terms go to except 
the last term which goes to q. Therefore, C/„ — >■ g as n — >■ cx). 



Observe that 

IKI<E 



x{DXr.Y \\Xn{DX„ 



< WE 



X (DXn) 



E 



\Xn{DXn)\\; 



= WE 



x2 {D.XnT df^ (z) 



E 



l^n(^^n)ll;, 



Note that 



E 



\Xn {DXn 



E 



^'ll^^nll', 



Un q- 



From the third point, E 
Finally for W„, 

|W„| =E 



— > so as n — )• 00. 



[D^Xn) d^i (z) 



< E 



x2 {D,Xnfdii{z) 



Putting them together we get the fourth point: E[X^] ^ 3 as n 00. ■ 

In Wiener space, convergence in a fixed Wiener chaos to a standard normal distribution is characterized 
by the convergence of the fourth moments to 3 or of the convergence of the norm of certain contractions to 
(see list preceding Corollary [^T]). We would then like to see if the same situation holds in Wiener-Poisson 
space. At this point, this appears to be an open question. We then finish with the following theorem which 
shows convergence in distribution and of the fourth moments to 3 if certain contractions converge to 0. 
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Theorem 29 Let Xn = Iq (fn), with E [X^] = ql |l/n||^»g — >■ 1- Assume that for r — 0, . . . ,q — 1 and 
s = 0,...,q-r, ||/„ (g)f, fn\\f^»(2g-2^-.) l{s=o,r#o}u{s#o,r=o} ^ 0. Then as oo, 

• Xn — > A/'(0,1) m distribution. 

Proof. The first assertion is the fourth point in Proposition [^S) For the second point, we refer to Corollary 

^ in L^{n) as 



BTlto see that it suffices to prove (X„) — gx„ — > in (fl) and (^\x\ {DX^f , \DX. 
n ^ oo. These are immediate when we note that 



9*{Xn) - gx„ ^ ^ ll^-'^nlll in ^^(ri) (by point 2 of Proposition 



and 



E 



\x\ {DXnf , \DX„ 



< E 



xiDXnr \\DX„ 



< WE 



X {DXnY 



E 



\DXr, 



E 



x2 {DXnf dn 



^/qE\X^ (by point 3 of Proposition 
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